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I  INTRODUCTION 


Under  the  sponsorship  of  the  Defense  Nuclear  Agency  (Contract 
DNA001-75-C-0206,  SRI  International  successfully  fielded  a  UHF/SHF 
coherent  transmission  experiment.  This  experiment  was  designed  to  measure 
the  effects  of  dust  and  debris  on  signals  passing  through  the  cloud 
lofted  by  the  DICE  THROW  high-explosive  test.  The  test,  which  was  a 
628-ton  ANFO  (ammonium  nitrate  plus  fuel  oil)  detonation,  took  place  at 
0800  MDT  on  6  October  1976  at  the  Giant  Patriot  site  in  the  northern  part 
of  the  White  Sands  Missile  Range  (WSMR)  (Figure  1).  Because  the  explosion 
lofted  between  4  and  10  kt  of  (desert)  soil,  it  simulated  in  at  least 
one  way  a  muvh  larger  nuclear  detonation  than  its  500-ton  TNT  equivalent. 

A  rather  complete  description  of  that  test  and  the  various  experimental 
programs  (including  this  one)  associated  with  it  appears  in  the  Proceedings 

ic 

of  the  DICE  THROW  Symposium.1 

The  UHF/SHF  transmission  experiment  consisted  basically  of  making 
measurements  of  amplitude  and  phase  fluctuations  on  12  different  CW  sig¬ 
nals  between  379  MHz  and  10,2  GHz.  Those  signals  were  transmitted  along 
six  paths  penetrating  the  space  above  and  in  the  vicinity  of  ground  zero 
where  the  dust  cloud  was  expected.  Preliminary  results  of  the  experiment 
and  a  report  on  investigations  of  soil  properties  appear  in  topical  re¬ 
ports  1  and  2  under  this  contract. 3* 15  An  earlier  report4  covered  the 
experimental  design  phase. 

The  primary  reason  for  fielding  this  experiment  and  making  these 
measurements  was  simply  that  until  then,  no  one  had  made  any  rf  trans¬ 
mission  measurements  through  high-density  dust  and  debris  clouds  lofted 
by  large  explosions.  There  is  legitimate  concern  that  various  kinds  of 
communication  and  sensor  systems  might  be  affected  adversely  by  large- 
scale  dust  clouds  produced  by  surface,  or  near-surface  nuclear  detonations. 


"k 

References  are  listed  at  the  end  of  the  report. 
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Our  goal  was  to  make  accurate  measurements  of  both  amplitude  and  phase 
perturbations  ov*r  a  range  of  wavelengths  that  are  commonly  used  by  com¬ 
munications  and  sensing  systems.  Although  many  of  the  frequencies  used 
were  in  or  near  bands  used  by  military  systems,  there  was  no  attempt  to 
simulate  in  any  way  the  operation  of  any  actual  system.  These  measure¬ 
ments  were  strictly  of  basic  rf  effects  on  simple,  CW,  signals.  Nuclear 
dust  cloud  rf  effects  can  in  principle  be  estimated,  using  these  basic 
measurements. 

Of  course,  the  magnitudes,  durations,  and  extents  of  the  rf  effects 
measured  during  DICE  THROW  were  much  smaller  than  are  likely  to  be  en¬ 
countered  in  a  nuclear  environment.  Although  the  DICE  THROW  main  event 
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was  much  more  effective  at  raising  dust  than  the  "megaton  of  dust  per 
megaton  of  yield"  rule  would  predict,  it  is  still  quite  small,  and  the 
flow-field  velocities  for  the  DICE  THROW  charge  configuration  and  composi¬ 
tion  are  very  small  as  well.  According  to  calculations,  the  capped  cylin¬ 
der  ANFO  configuration  produces  flow  fields  having  significantly  lower 
velocities  than  even  those  occurring  after  tangent-sphere  TNT  detonations. 
In  addition,  the  areal  extent  of  the  dust  cloud  was  only  on  the  order  of 
a  few  hundred  meters  during  the  early  phases.  After  a  nuclear  detonation, 
a  much  larger  area  and  volume  are  affected,  and  the  larger  velocities 
during  the  rising  fireball  phase  would  tend  to  support  more  and  larger 
particles  for  much  greater  lengths  of  time.  Extrapolation  to  the  nuclear 
environment,  while  not  part  of  the  work  reported  here,  is  an  important 
and  in  many  ways  a  difficult,  uncertain  step  that  needs  to  be  taken. 

These  data  provide  one  of  the  bases  for  such  extrapolation. 

The  first  subject  in  this  report  if.  a  description  of  the  experiment 
and  its  layout.  The  next  two  sections  are  devoted  respectively  to  pre¬ 
senting  and  to  interpreting  the  measured  effects,  which  can  be  divided 
into  two  categories.  The  first  category,  which  lasted  for  several  seconds 
following  the  detonation,  was  characterized  by  intense  fading,  including 
the  complete  loss  of  a  coherent  signal  on  many  of  the  frequencies.  Scat¬ 
tering  from  large  particles  seems  to  have  been  very  important  then.  The 
second  category  was  characterized  by  gentle  effects,  because  of  fallout 
of  the  larger  particles.  The  duration  of  effects  was  limited  by  the  rapid 
drift  of  the  dust  cloud  away  from  the  signal  paths  due  to  higher-than- 
expected  winds  at  test  time. 

There  are  two  appendices.  The  first  (Appendix  A)  is  a  brief  review 
ci  the  theoretical  effects  of  dust  clouds  on  electromagnetic  propagation 
in'  the  JHF/SHF  range.  Results  of  laboratory  measurements  of  the  dielec¬ 
tric  properties  of  soil  samples  gathered  from  the  DICE  THROW  crater  are 
summarized  there  as  well.  Our  conclusion  is  that  phase  shifts  and  attenu¬ 
ation  caused  by  absorption  are  probably  the  dominant  effects  except  at 
early  times  before  sufficient  numbers  of  large  particles  have  fallen  out. 
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Scattering  losses  can  become  very  large  when  very  many  large  particles 
are  in  the  dust  cloud.  These  results  are  used  in  the  interpretation  o£ 
the  experimental  results. 

Development  of  the  dust  cloud  is  considered  in  Appendix  B,  Photo¬ 
graphic  records,  some  of  which  were  obtained  from  our  own  boresight 
cameras,  are  compared  with  other  events  and  with  Air  Force  Weapons  Labora¬ 
tory  (AFWL)  hydrodynamic  (HULL)  code  calculations.  Comparisons  are  also 
made  with  some  previous  KE  test  results,  and  with  the  appearance  of  several 
nuclear  events.  These  comparisons  illustrate  the  substantial  dissimilari¬ 
ties  between  the  clouds  produced  by  different  types  of  detonations. 

The  measured  DICE  THROW  crater  volume  was  used  to  bound  the  uncertainty 
in  the  amount  of  lofted  dust  and  debris.  A  simple  model  for  the  average 
dust  density  that  is  based  on  the  measured  dust  cloud  volume  and  on  esti¬ 
mates  of  the  mass  of  lofted  soil  is  developed  from  the  photographs. 

These  photographic  records  are  very  important  to  the  interpretation  of 
the  results  of  the  experiment. 


II  EXPERIMENTAL  LAYOUT  AND  EQUIPMENT 


Figure  2  is  a  schematic  of  the  experimental  layout.  Signals  were 
transmitted  from  two  locations  and  were  received  at  three  sites,  provid¬ 
ing  six  measurement  paths  (labeled  Paths  1  through  6).  In  addition,  the 
413-MHz  signal  from  the  main  transmitter  was  received  at  the  "phase- 
repeater  station,"  multiplied  to  1239  MHz,  then  retransmitted  to  the 
receiving  system.  This  arrangement  served  to  transfer  a  phase  reference 
signal  from  the  main  transmitter  to  the  receiving  system  along  a  path 
unaffected --or  only  briefly  affected--by  the  explosion.  Finally,  a 
1274-MHz  signal  from  the  main  transmitter  was  received  at  North  Oscuro 
Peak  to  phase-lock  the  remote  transmitter.  The  transmitter  power  levels 
ranged  between  10  and  100  mW. 

Table  1  lists  the  exact  frequencies  of  the  signals  transmitted  along 
the  six  measurement  paths.  Initially,  no  measurements  were  planned  for 
Paths  5  and  6.  It  was  possible,  however,  to  receive  the  424-MHz  signal 
from  the  remote  transmitter  at  any  of  the  receiving  antenna  sites.  There 
fore,  because  it  appeared  that  a  net  gain  of  data  could  easily  be  ob¬ 
tained,  we  decided  to  make  a  simple  modification  to  switch  sequentially, 
or  commutate,  the  424-MHz  receiver  channel  between  the  three  antenna 
sites.  A  2-s  period  was  chosen,  with  the  424-MHz  UHF  signal  of  Path  4 
received  for  1  s  and  the  424-MHz  UHF  signals  of  Paths  5  and  6  received 
for  0.5  s  each.  Although  Path  6  passed  quite  far  from  ground  zero,  it 
was  believed  that  there  was  a  reasonable  chance  that  the  cloud  might 
drift  into  that  path  later. 

Two  major  considerations  primarily  dictated  the  physical  locations 
of  the  various  components  of  the  experiment.  First,  based  on  the  ex¬ 
pected  size  of  the  dust  cloud,  a  20-m  Fresnel  zone  radius  at  413  MHz  was 
desirable.  Second,  Path  4  was  required  to  pass  about  200  m  directly  over 
ground  zero.  Secondary  considerations  included:  (1)  access  and  practica 
bility,  (2)  keeping  the  signals  clear  of  obstacles  and  equipment,  and 
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REMOTE 


FIGURE  2  EXPERIMENTAL  LAYOUT  OF  OICE  THROW  TRANSMISSION  EXPERIMENT 


Table  1 

MEASUREMENT  FREQUENCIES 


(3)  separating  Paths  2  and  3,  respectively,  by  one  (40  m)  and  three 
(120  m)  UHF  Fresnel  zone  diameters  from  ground  zero  (GZ)  in  order  to 
probe  the  dust  cloud  laterally  along  independent  lines  of  sight  at  UHF. 
Figures  3  and  4,  respectively,  present  the  plan  and  elevation  (in  the 
plane  of  ground  zero)  views  of  the  resulting  layout.  The  WSMR  coordi¬ 
nates  of  the  various  sites  are  listed  in  Table  2.  (So  that  it  would  be 
more  accessible,  the  phase  repeater  station  was  not  located  at  the  sur¬ 
veyed  point,  and  because  its  new  position  was  not  resurveyed,  only  its 
approximate  location  in  WSMR  coordinates  is  listed.)  For  reference,  the 
WSMR  coordinates  of  ground  zero,  in  meters,  were  E135,027,  N198,205,  and 
H1442  (see  Figure  3), 

The  beamwidths  of  the  antennas  were  chosen  so  that  the  entire  ex¬ 
pected  extent  of  the  early  dust  cloud  would  fall  well  within  the  beams. 
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ORTH  OSCURO 
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FIGURE  3  DICE  THROW  SITE  LAYOUT  WITH  MICROWAVE  TRANSMISSION  EXPERIMENT  SITE  LOCATIONS 


PATH  4^ 


PATH  5 


PATH  6 


PATH  1 


LPATH  2  f 
r  10-11  m 


PATH  3 


LA-3972-6 

FIGURE  4  LINE-OF-SIGHT  PENETRATION  POINTS  IN  VERTICAL  PLANE  THROUGH 
GROUND  ZERO  AND  PERPENDICULAR  TO  PATHS  1  AND  4 


This  choice  tended  to  eliminate  any  antenna  beam-limiting  effects  on  the 
measurements.  However,  the  antenna  beams  had  to  be  as  narrow  as  practical 
in  order  to  reduce  the  possibility  of  troublesome  reflections  from  nearby 
objects.  Except  for  UHF,  where  impractically  large  antennas  would  have 
had  to  be  used,  the  antenna  parameters  as  listed  in  Table  3  fulfilled 
these  experimental  requirements,  An  extensive  program  of  antenna  measure¬ 
ments  was  undertaken  before  the  experiment  was  fielded. 

Because  the  main,  multi frequency  transmitter  had  to  be  placed  quite 
close  (630  m)  to  ground  zero  in  order  to  reduce  the  UHF  Fresnel  zone,  it 
was  necessary  to  harden  the  antennas  and  to  protect  the  electronics  from 
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Table  2 


WSMR  COORDINATES  OF  DIGE  THROW  MICROWAVE 
TRANSMISSION  EXPERIMENT  SITE  LOCATIONS 


WSMR  Coordinates 
(m) 

Site 

E 

N 

H 

Main  Transmitter 

135,369 

198,525 

1450 

Remote  Transmitter 
(North  Oscuro  Peak) 

148,810 

206,340 

2440 

Phase  Repeater 

135,600 

196,200 

1440 

Receiving  Site  1 
(Main) 

131,614 

196,191 

1429 

Receiving  Site  2 
(Outrigger  1) 

131,843 

195,995 

1429 

Receiving  Site  3 
(Outrigger  2) 

131,972 

195,319 

1428 

the  shock  wave.  (The  measured  overpressure  at  630  m  was  about  19  kPa.) 
Protection  and  immobilization  of  the  antennas  was  accomplished  by  casting 
them  into  a  concrete  block  and  filling  their  apertures  with  foam.  The 
UHF  antenna,  which  had  a  wider  beamwidth  than  the  others,  was  skewed  to 
the  left  in  order  to  increase  the  signal  strength  at  the  phase  repeater 
station  (the  others  were  aimed  in  azimuth  at  ground  zero).  An  electronics 
shelter  was  partly  buried  in  the  rear  of  a  20-ft  earthen  mound  constructed 
to  support  and  elevate  the  antennas.  No  damage  was  sustained  by  either 
the  antennas  or  the  electronics  during  the  evenc.  The  phase  repeater 
station  consisted  of  antennas,  electronics,  and  cameras  mounted  on  a 
20-ft  tower. 

Several  actions  were  taken  to  reduce  the  possibility,  and  the  ef¬ 
fects,  of  any  multipath  signals  due  to  reflections  from  the  ground  surface 
or  from  equipment  or  structures.  Circular  polarization  was  used  on  all 
of  the  measurement  paths  because  specular  or  quasi-specular  reflection 
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from  conductive  surfaces  reverses  the  sense  of  circular  polarization. 
(Because  the  requirements  were  much  less  severe  on  the  phase-reference 
paths,  vertical  linear  polarization  was  used  on  them.)  Both  the  trans¬ 
mitting  and  receiving  antennas  were  tilted  upward  by  one-half  of  their 
Ibeamwidths  to  discriminate  against  ground -reflected  signals.  Clutter 
fences  were  placed  in  front  of  the  receiving  antennas  to  reduce  ground 
reflections.  These  fences  were  30  m  wide  and  were  placed  at  the  geo¬ 
metrical  points  of  reflection  (50  m  in  front  of  the  receiving  antennas) 
for  surface-reflected  signals  from  the  North  Oscuro  transmitter.  Although 
such  fences  are  quite  effective  against  surface-reflected  signals  from 
high-elevation  sources,  they  are  not  very  effective  in  suppressing  those 
associated  with  very-low-elevation-angle  paths  like  Paths  1  to  3  because 
the  diffraction  angles  are  very  small.  The  best  one  could  do  would  be 
to  place  the  top  of  the  fence  so  that  both  the  direct  and  surface- 
reflected  signals  would  be  in  the  shadow  zone;  such  an  arrangement  at¬ 
tenuates  any  reflected  signals  more  than  it  does  the  direct  one.  For 
DICE  T'lROW  the  edges  of  the  geometric  shadows  of  the  tops  of  the  2.3-ro- 
high  (7.5  ft)  fences  illuminated  by  the  main  transmitter  passed  through 
the  centers  of  all  of  the  receiving  antennas  behind  the  fences.  All  of 
the  receiving  antennas  were  aimed  at  the  appropriate  transmitters.  The 
direct  signals  from  the  upper  paths  (4  through  6)  were  virtually  unaf¬ 
fected  by  the  clutter  fences. 

The  receiver  was  designed  to  be  linear  over  at  least  a  60-dB  dynamic 
range,  and  the  signal  levels  were  adjusted  to  be  about  10  dB  below  satura¬ 
tion  to  account  for  any  enhancements  due  to  diffraction  or  focusing. 

Rather  than  detecting  amplitude  and  phase,  the  quadrature  signal  compo¬ 
nents  (which  can  be  interpreted  as  the  real  and  imaginary  parts  of  the 
complex  representation  of  the  received  signal  or  as  a  pair  of  zero- 
carrier-frequency  IF  signals)  were  developed  at  the  receiver  outputs. 

The  amplitude  and  phase  of  each  signal  can  be  easily  computed  from  each 
pair  of  quadrature  components.  In  a  sense,  the  nonlinear  "detection" 
process  occurs  in  the  mathematical  computation.  Greater  processing 
flexibility  justifies  the  quadrature-component  approach.  The  overall 
bandwidths  of  a] 1  the  receiver  channels  were  limited  at  each  quadrature 
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component  output  by  a  four-pole,  low-pass  Butterworth  filter  with  150-Hz,. 
3-dB  cutoff  frequency.  Each  of  the  13  receiver  channels  used  provided 
a  pair  of  quadrature  component  outputs.  One  of  the  13  channels  was  used 
for  the  1239-MHz  phase-reference  signal.  One  of  the  twelve  recording 
channels  used  for  measurement  paths  was  commutated  between  the  424-MHz 
signals  from  Paths  4,  5,  and  6  as  described  above.  Quadrature  component 
pairs  were  simultaneously  digitized  by  a  pair  of  12-bit  (±2048-point) 
analog-to-digital  converters  (one  for  each  of  the  quadrature  components). 
Thus,  26  channels  of  data  were  digitized.  The  digitization  rate  was  500 
samples  per  second,  and  the  commutation  interval  was  80  us.  Then,  under 
the  control  of  an  HP- 2100  minicomputer,  the  data  were  written  on  magnetic 
tape.  Two  independent  digital  tape  recorders  were  used.  A  backup  analog 
recording  system  was  also  operated.  That  system  multiplexed  a  50-kHz 
IF  signal  from  each  of  the  13  receiver  channels  with  a  100-kHz  pilot  tone 
serving  as  a  phase  reference,  and  recorded  each  resulting  composite  sig¬ 
nal  on  a  single  tape  recorder  track.  Time  signals  were  also  recorded  by 
both  the  primary  and  the  backup  recording  systems, 
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Ill  MEASURED  RF  EFFECTS 


A.  General 

In  this  section  we  present  the  final  results  of  the  rf  transmission 
experiment.  The  data,  which  were  originally  presented  in  Ref.  1,  have 
been  carefully  calibrated  and  purged  of  all  the  artifacts  that  could  be 
identified.  A  particularly  great  effort  was  expended  toward  the  goal  of 
obtaining  continuous  phase  records--i.e. ,  removing  the  discontinuities 
at  ±180°.  While  this  waf  a  simple  process  conceptually,  in  practice  it 
was  very  difficult  to  devise  an  algorithm  to  make  the  decision  of  whether 
or  not  a  phase  transition  across  the  ±180°  boundary  actually  occurred 
when  the  signals  were  weak  and  fading. 

The  duration  of  rf  effects  was  quite  brief.  On  the  lower  paths 
(1,  2,  and  3),  there  was  a  period  of  very  strong  effects  lasting  until 
about  T  +  3toT+4s.  ByT  +  25s,  all  signal  perturbations  had  ceased 
on  those  paths.  Of  the  upper  paths,  only  Path  4  passed  through  the  dust 
cloud  because  of  the  adverse  winds,  (Although  we  did  not  expect  an  occulta- 
tion  on  Path  6  due  to  its  distance  from  GZ,  we  believed  that  there  was  a 
good  chance  of  obtaining  one  on  path  5.)  Rf  perturbations  on  Path  4 
lasted  until  T  +  37  s,  when  the  dust  cloud  blew  out  of  the  line  of  sight. 

For  a  quick  comparative  look  at  the  evolutional  growth  of  the  DICE 
THROW  cloud  as  seen  from  the  three  major  camera  sites,  R-l,  PR,  and  NOP, 
Figures  5  through  10  have  been  prepared.  Figure  5  displays  the  cloud 
outlines  as  recorded  on  film  from  the  R-l  site  at  0.1,  4,  10,  and  20  s 
superimposed  on  a  photograph  taken  at  T  +  0.1  s,  which  is  lined  up  with 
its  corresponding  outline.  In  addition,  the  propagation  path  intersections 
are  overlaid  in  their  correct  locations.  Figure  6  similarly  displays  R-l 
contours  at  20,  30,  60,  and  90  s  over  a  photograph  at  20  s. 

Figures  7  and  8  show  contours  overlaid  on  the  earliest  photograph 
for  each  given  time  sequence  as  well  as  raypath  configuration  as  seen 


PATH  INTERSECTION  POINTS 


PATH  INTERSECTION  POINTS 


EVOLUTIONAL  GROWTH  OF  OICE  THROW  CLOUD  FROM  PR  SITE  AT  T  +  0.2,  4,  10,  AND  20  $  WITH 
PROPAGATION  PATHWAYS 


FIGURE  8  EVOLUTIONAL  GROWTH  OF  DICE  THROW  CLOUD  FROM  PR  SITE  AT  T  +  20,  30,  AND  60  s  WITH  PROPAGATION 
PATHWAYS 


PATH  INTERSECTION  POINTS 


FIGURE  10  EVOLUTIONAL  GROWTH  OF  DICE  THROW  CLOUD  FROM  NOP  SITI 
AND  3  min  WITH  PROPAGATION  PATH  INTERSECTION  POINTS 
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from  the  FR  site.  Figures  9  and  10  present  the  view  of  the  other  side 
of  the  cloud  seen  from  NOP.  In  the  second  figure  for  each  site  of  the 
series,  the  early-time  contour  is  the  latest  time  contour  outlined  in 
the  first  figure. 


B.  Time  Histories  of  RF  Effects 

Figures  11  and  12  show  the  amplitudes  and  phases  for  the  lower  path 
UHF  signals  from  T  -  1  to  T  +  14  s.  The  sign  of  the  phase  shift  produced 
by  the  dust  cloud  for  each  of  the  signals  depends  on  the  manner  in  which 
the  signal  is  translated  in  frequency  by  the  receiver,  and  can  be  either 
positive  or  negative.  For  all  of  these  UHF  signals,  the  dust  cloud  would 
produce  negative  phase  shifts.  The  Path  1  phase  shifts  of  Figure  11  indi¬ 
cate  that  there  was  a  relatively  uniform  fallout  of  most  of  the  material 
in  the  cloud  that  ended  at  about  T  +  5  s.  No  particular  change  in  the 
appearance  of  the  cloud  accompanied  the  end  of  the  rapid  fallout  phase. 

It  should  be  noted  that  a  return  to  the  baseline  signal  level  could  not 
be  expected  in  general  because  of  the  blocking  effect  of  the  explosive 
stack. 
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It  is  possible  that  the  passage  of  the  shock  wave  at  various  places 
could  perturb  the  signals  in  undesirable  ways,  such  as  introducing 
transient  and  permanent  phase  shifts  with  amplitude  changes.  Table  4 
lists  the  shock-wave  arrival  times  at  several  locations.  These  times 
were  obtained  from  the  average  of  the  measured  arrival  times  along  the 
three  instrumentation  gauge  lines  laid  out  for  the  experiment.1  The 
antenna  oscillations  at  the  main  receiver  site  and  the  direct  effects  of 
the  shock  waves  passing  through  the  signal  paths  are  the  only  ones  that 
have  been  identified. 

Because  the  rf  effects  of  the  dust  cloud  were  very  much  the  same  on 
the  Path-1  UHF  signals,  which  were  spaced  69  MHz  apart  (Figure  11),  the 
conerent  bandwidth  at  400  MHz  was  substantially  wider  than  70  MHz.  Since 
the  coherent  bandwidth  is  inversely  related  to  the  multipath  time  spread, 
this  result  is  consistent  with  the  small  lateral  extent  of  the  dust  cloud. 
Strong  diffractive  effects,  however,  were  clearly  present.  Distortion  of 
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FIGURE  12  UHF  SIGNAL  FLUCTUATIONS  FOR  PATHS  1,  2,  AND  3 
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Table  4 

SHOCK-WAVE  ARRIVAL  TIMES  AT  VARIOUS  LOCATIONS 
B.  SED  ON  MEASURED  FREE -FIELD  BLAST  DATA 


Position 

Range 

<m) 

Arrival  Time 

(sj 

Main  Transmitter  (MT) 

630 

1.15 

Phase  Repeater  (PR) 

2400 

0.8 

Receiver  (R-l) 

3900 

k 

11.0 

Perpendicular  to 

MT-PR  path 

535 

0.9 

Perpendicular  to 

PR  to  R-l  path 

2330 

6.5 

* 

Extrapolated, 


the  phase  shifts  due  to  diffraction  was  severe  enough  to  preclude  inter¬ 
pretation  of  phase  shift  in  terms  of  integrated  dust  density  until  times 
later  than  T  +  1  to  T  +  1.5  s. 

This  type  of  distortion  is  very  evident  in  the  results  of  Path  2 
between  Tq  and  T  +  1.5  s  shown  in  the  center  panel  of  Figure  12.  Part 
of  the  positive  phase  shift,  which  implies  a  loss  of  matter,  in  that  time 
interval  could  have  been  caused  by  a  decrease  in  atmospheric  refractivity 
due  to  heating.  As  will  be  seen,  this  effect  did  occur  along  Path  4.  If 
it  were  significant  here,  it  would  tend  to  reduce  the  effects  of  the  dust. 
For  instance,  the  HULL  calculations  predict  temperatures  above  600  K 
within  a  150-m  radius  around  ground  zero  near  the  surface  during  the 
first  three  seconds  after  detonation.  Since  the  dominant  term  in  the 
expression  for  atmospheric  refractivity  is  inversely  dependent  on  tempera¬ 
ture,  these  temperatures  would  reduce  the  refractivity  on  the  order  of 
100  N-units.  At  400  MHz,  the  corresponding  phase  shift  over  a  300-m  path 
is  about  15°.  Because  this  value  is  less  than  10%  of  the  measured  phase 
shifts  at  early  times  on  Path  1,  temperature  effects  were  probably  not 
very  important,  although  irregularities  may  have  contributed  to  some  of 
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the  fluctuations  observed.  The  distortion  due  to  the  diffraction  on 
Path  2  was  so  great  that  the  net  phase  shifts  were  much  larger  -in  the 
positive  direction  than  can  be  explained  by  heating. 

The  very  strong  similarities  between  the  amplitude  pattern  for  Path  2 

from  T  +  0.5  s  to  T  +  1.5  s  and  those  for  Path  1  from  T  to  T  +  1  s  lead 

o 

to  the  conclusion  that  the  Path  1  phase  shifts  were  distorted  in  the  same 
way  as  that  for  Path  2.  This  means  that  the  undistorted  Path  1  phase 
shifts  were  probably  larger  than  those  measured. 

Another  propagation  effect  is  that  a  dust  cloud  tends  to  act  as  a 
converging  lens.  This  is  the  reason  for  the  early-time  10-dB  signal 
enhancement  on  Paths  1  and  2.  Note  that,  although  the  visible  dust  cloud 
practically  instantaneously  enveloped  all  of  the  lower  paths,  the  focus 
for  Path  2  was  delayed  about  0.5  s  behind  those  for  Path  1.  Therefore 
the  focus  is  apparently  due  to  a  dense  bulk  of  crater  ejecta  within  the 
visible  cloud.  Since  the  paths  were  spaced  40  m  apart,  the  0.5-s  delay 
implies  a  horizontal  velocity  of  about  80  m/s  for  the  dense  ejecta  bulk. 
Although  significant  extinction  would  be  expected  from  a  very  dense  region, 
phase  effects  were  dominant  and  the  rf  energy  at  400  MHz  was  focused  much 
more  strongly  than  it  was  attenuated. 

As  would  be  expected,  the  onset  of  the  strongest  effects  on  Path  3 
were  delayed  even  more.  The  top  panel  of  Figure  12  shows  that  the  in¬ 
tensity  of  the  fluctuations  increased  at  T  +  2  s  and  changed  somewhat  in 
character.  However,  for  this  path  there  was  no  general  phase  shift  due 
to  the  integrated  dust  density.  Rather,  the  signal  appears  to  have 
exhibited  Rician  fading.  Figure  13(a),  which  is  a  plot  of  the  quadrature 
components  of  the  413-MHz  Path-3  signal  between  Tq  and  T  +  5  s,  supports 
this  interpretation.  The  signal  seems  to  be  rather  well  characterized  as 
being  the  sum  of  a  steady,  "coherent"  component  and  a  random,  "scattered" 
component  having  no  preferred  phase;  in  contrast,  the  Path  1  signal, 
whose  quadrature  components  from  Tq  to  T  +  5  s  are  shown  in  Figure  13(b), 
had  a  coherent  component  that  underwent  a  large  phase  change.  This  means 
that  the  part  of  the  dust  cloud  penetrated  by  Path  3  was  very  inhomogeneous 
and  did  not  have  a  very  high  average  density.  The  2-s  delay  and  120-m 
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offset  of  Path  3  from  ground  zero  implies  a  60-ra/s  ejecta  velocity.  This 
is  slightly  smaller  than  the  80-m/s  horizontal  velocity  deduced  from  the 
Path  2  delay,  but  is  consistent  with  a  lower  horizontal  velocity  for  the 
higher-angle  ejecta  that  would  occult  Path  3. 

The  higher-frequency  (L-,  S-,  and  X-band)  Path  1  signals  were  dif¬ 
ferent  from  the  UHF  signals  during  the  early  period  in  that  all  of  them 
exhibited  an  effectively  complete  loss  of  a  coherent  signal  component. 

Only  a  rapidly  fading,  scattered  component  of  the  transmitted  signal  was 
received  for  several  seconds  after  detonation,  during  which  time  phase 
data  were  nearly  meaningless.  Figure  14  shows  the  L-,  S-,  and  X-band 
amplitudes  and  phases  from  T  -  1  to  T  +  14  s.  Because  these  data  have 
been  smoothed  using  a  20-ms  running  average,  the  depth  of  the  fades  has 
been  suppressed  by  about  3  to  5  dB--i.e. ,  the  fades  were  about  3  to  5  dB 
greater  in  the  raw  data.  The  loss-of-coherent-signal  periods  manifest 
themselves  in  Figure  14  in  three  ways:  (1)  a  decrease  in  the  average 
signal  level,  (2)  rapid  and  deep  fading,  and  (3)  a  tendency  for  the  phase 
to  wind  up  by  many  revolutions,  followed  by  a  recovery  to  a  new  baseline. 
There  is  no  explanation  for  the  tendency  for  the  phase  to  wind  up  in  the 
negative  direction;  it  may  be  an  artifact.  A  change  in  the  baseline  value 
was  expected  because  of  the  disappearance  of  the  charge,  which  was  in  the 
Path  1  line  of  sight. 
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Because  dust  cloud  effects  become  stronger  as  the  frequency  increases, 
recovery  of  a  sufficiently  strong  coherent  component  to  make  interpreta¬ 
tion  of  phase  shifts  in  terms  of  integrated  dust  density  possibly  takes 
longer  at  higher  frequencies.  The  Path  1  L-band  signal  recovered  at  about 
T  +  2  s,  while  the  S-  and  X-band  recoveries  were  delayed  until  T  +  2.5  s 
and  T  +  3  s,  respectively.  After  recovery,  the  phases  relaxed  to  their 
new  baseline  values  in  a  way  consistent  with  the  VHF  results. 

At  about  T  +  11  s  the  shock  wave  arrived  at  the  receiver  and  per¬ 
turbed  the  system  enough  to  knock  the  X-band  local  oscillator  out  of 
phase  lock.  The  noise  level  was  10  to  13  dB  below  the  signal  level  during 
the  intense  fading  period  Tq  to  T  +  3  s,  which  proves  that  although  the 
coherent  signal  disappeared,  there  was  some  signal  getting  through.  Note 
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FIGURE  14  L-,  S-,  AND  X-BAND  AMPLITUDES  AND  PHASES 

FROM  T  -  1  s  TO  T  +  14  s 
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that  there  was  no  tendency  for  the  phase  to  wind  up  in  any  particular 
direction  when  there  was  only  noise  preset.  The  tendency  for  the  phase 
to  wind  up  in  a  preferred  direction  during  the  intense  fading  period  was, 
therefore,  probably  a  true  dust-cloud  effect  rather  than  an  instrumental 
artifact. 

Figure  15  is  a  polar  plot  of  the  S-  and  X-band  quadrature  components 
from  Tq  to  T  +  5  s.  At  S-band,  before  the  recovery  phase,  which  was 
characterized  by  a  spiral  trajectory,  the  signal  appears  to  be  distributed 
in  a  manner  consistent  with  Rayleigh  statistics.  The  L-  and  X-band  signals 
behaved  in  similar  ways  except  that  the  spiral  recovery  trajectory  in  the 
X-band  case  (Figure  15)  cannot  be  distinguished  because  it  performed  seve¬ 
ral  revolution:!.  Statistical  computation  also  seem  to  support  Rayleigh 
statistics  at  early  times.  Table  5  presents  the  results  from  computing 
various  moments  from  the  data  in  several  time  intervals.  During  the 
0-to-2.7-s  interval,  the  mean  values  were  much  smaller  than  the  standard 
deviations,  which,  furthermore,  were  approximately  equal  to  one  anocher. 
Also,  the  correlation  coefficient,  r,  between  I  and  Q  was  very  small. 
However,  these  results  may  be  very  misleading,  for  we  know  that  the  phase 
behaved  in  a  very  deterministic  way  during  that  interval.  Table  6  lists 
the  ratio  of  the  variance  of  the  L-,  S-,  and  X-band  signal  power  to  its 
mean  value  for  the  0-to-2.7-s  interval.  If  the  signals  did  follow  Ray¬ 
leigh  statistics,  this  ratio  would  have  been  unity.  Because  the  actual 
S-  and  X-band  ratios  were  very  different  from  that  value,  we  must  conclude 
that  those  signals  were  not  fading  in  a  Rayleigh  fashion.  In  fact,  they 
exhibited  much  greater  variability  than  expected  for  the  Rayleigh  case, 
particularly  at  S-band. 

The  Path  2  and  Path  3  X-band  signals  also  experienced  periods  of 
intense  fading,  during  which  phase  changes  cannot  be  related  directly  to 
integrated  dust  density.  Consistent  with  the  progressive  delay  of  the 
strongest  effects  with  path  offset  from  ground  zero  seen  in  the  UHF  dnta 
(Figure  12),  the  Path  2  and  3  X-band  data  are  also  progressively  delayed 
(Figure  16).  The  strongest  effects  on  Path  2  spanned  the  interval 
T  +  1.5  to  T  +4.5  a,  and  those  on  Path  3  occurred  between  T  +  4.5  s  and 
T  +  8  s.  This  pattern  implies  that  the  intense  fallout  of  larger 
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Table  5 

PATH  1  QUADRATURE  SIGNAL  COMPONENT  STATISTICS 
IN  UNITS  OF  DIGITIZER  OUTPUT  COUNTS 


Time  Interval 
(s) 

i 

Q 

°I 

°Q 

0-2.7 

65 

67 

EH| 

480 

2. 7-5. 5 

-150 

-1100 

EM 

310 

5-8 

960 

-980 

310 

320 

0-2.7 

-49 

31 

250 

240 

2. 7-5.5 

290 

-490 

980 

680 

5-8 

880 

920 

700 

760 

0-2.7 

-5.4 

-5.1 

70 

68 

2. 7-5. 5 

| 

32 

82 

260 

270 

5-8 

260 

2.7 

480 

550 

0.005 

0.066 

-0.83 

0.004 

0.14 

0.29 


Table  6 

RATIOS  OF  SIGNAL  POWER  VARIANCE  (a  ) 

TO  MEAN  SIGNAL  POWER  (P) 

IN  THE  T  -to-2. 7-s  TIME  INTERVAL 
o 


Band 

p 

°P 

Op/P 

L 

5.0  X  105 

4.9  x  105 

0.98 

S 

1.2  X  105 

3,6  x  105 

3.00 

X 

9.5  X  103 

1.4  x  104 

1.45 

particles  within  the  dust  cloud  struck  the  ground  in  roughly  the  form  of 
a  slowly  expanding  annulus  about  80  to  95  m  wide  and  moving  outward  from 
ground  zero  at  about  27  m/s.  The  statistical  nature  of  these  signals 
during  their  respective  intense  fading  periods  paralleled  that  of  the 
Path  1  signal,  and  leads  to  the  same  conclusion  that  although  the  fadings 
were  similar  in  some  ways  to  the  Rayleigh  case,  they  were  much  more 


variable.  This  may  be  due  to  temporal  nonstationarity  of  conditions 
within  the  dust  cloud  and  changes  in  the  path  lengths  through  the  cloud 
with  time. 

Consistent  with  the  cloud  motion  (Figures  5,  6,  9,  and  10),  the  rf 
effects  on  Path  1  ceased  at  T  +  25.5  s.  Figure  17  shows  the  Path  1 
S-band  amplitude  and  phase  between  T  +  3  and  T  +  30  s.  Following  the 
initial  fallout,  the  phase  shift  was  relatively  steady  at  about  -35°  from 
T  +  12  to  T  +  21  s,  followed  by  a  fairly  linear  return  to  the  clear - 
atmosphere  baseline  at  T  +  25.5  s.  Amplitude  fluctuations  were  rather 
weak  throughout  this  entire  period.  The  Path  2  and  3  behaviors  were 
similar  to  this  except  that  effects  terminated  at  T  +  22  and  T  +  15.5  s, 
respectively.  These  times,  plus  the  termination  time  of  Path  1,  corre¬ 
spond  closely  to  a  12.5-m/s  lateral  cloud  drift  speed.  The  Path  2  and  3 
* 

UHF  phase  shifts  were  about  -5  to  -6  before  starting  their  return  to 

their  clear-atmosphere  values.  These  values  are  consistent  with  the 

*=*-35°  value  measured  on  Path  1,  which,  using  the  average  of  the  values  in 

2 

Table  A-3,  corresponds  to  an  integrated  dust  density  of  about  2.7  g/cm  . 

The  strong  surface  wind  and  lower-than-anticipated  cloud  height  pre¬ 
cluded  all  but  a  grazing  occultation  on  the  upper  paths.  Figures  5  through 
10  show  that  only  a  small  part  of  the  dust  cloud  intercepted  Path  4. 

Paths  5  and  6  always  remained  clear  of  dust.  Figure  18  presents  the  8.9- 
GHz  Path  4  results.  For  this  signal,  an  increase  of  matter  in  the  path 
would  cause  a  positive  phase  shift;  therefore,  negative  phase  shifts  are 
due  to  the  heating  of  the  air.  The  long-term  background  phase  drift  is 
shown  as  a  dashed  line  in  Figure  18. 

According  to  the  photographic  records,  the  dust  cloud  did  not  inter¬ 
cept  Path  4  until  T  +  14.2  s.  All  of  the  effects  before  then  were  caused 
by  energy  deposited  by  the  shock  wave  and  by  the  heated  air  in  the  flow 
field  outside  of  the  dust  cloud.  This  is  consistent  with  the  HULL-code 
calculations,  which  predicted  that  the  dust  cloud  would  occupy  only  a 


The  X-band  signals  were  not  reacquired  until  about  T  +  58  s,  well  after 
any  effects  were  over. 
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FIGURE  18  PATH  4  X-BAND  SIGNAL  FLUCTUATIONS 


fraction  of  the  appreciable  flow  field  for  a  substantial  length  of  time 

(cf..  Figures  B-10,  B-15,  B-20,  and  B-24).  The  strong  phase  fluctuations 

in  the  interval  Tq  to  T  +  6  s  were  probably  due  to  turbulence  following 

the  shock  wave.  That  was  followed  by  a  period  lasting  more  or  less  until 

T  +  12  s,  when  the  phase  had  a  small  negative  value  resulting  from  heated 

air  in  the  signal  path  (the  oscillations  beginning  just  before  T  +  11  s 

were  caused  by  the  shock  wave  striking  the  receiving  antenna  structure). 

Either  the  volume  or  the  temperature  of  the  heated  air  increased,  starting 

around  T  +  12  s,  just  preceding  the  dust  cloud.  After  the  dust  cloud 

intercepted  the  signal  path,  beginning  at  T  +  14.2  s,  the  phase  began  to 

shift  in  a  positive  direction.  Both  dust  and  heating  effects  were  present 

during  the  occultation,  and  we  believe  that  there  would  have  been  close  to 

twice  the  approximately  60°  peak  positive  phase  shift  measured,  had  there 

been  no  heating  effects.  That  is,  the  phase  shift  due  to  dust  of  about 

+120°  was  partially  offset  by  a  heating  phase  shift  of  about  -60°.  This 

is  a  crude  estimate,  however,  because  there  is  no  direct  way  to  separate 

the  two  effects  and  it  is  based  on  the  pre-occultation  behavior.  In 

agreement  with  the  photographic  data,  the  phase  returned  to  its  background 

level  when  the  dust  cloud  cleared  Path  4  at  T  +  35  s.  For  phase  shifts 

of  60  and  120°  at  8.9  GHz  the  corresponding  integrated  dust  densities  are 
2 

1.75  and  3.5  g/cm  ,  respectively. 

Figure  19  shows  the  amplitude  and  phase  of  the  Path  4  UHF  signal, 
which  was  time-multiplexed  wit'.,  the  Path  5  and  6  signals.  Except  for  the 
opposite  sign  of  phase  changes,  the  effects  were  consistent  with  the 
S-band  results.  The  heating  effects  and  turbulence  caused  by  the  shock 
wave  perturbed  all  of  the  upper-path  signals  (Figure  20),  even  though 
Path  6  passed  720  m  away  from  ground  zero.  At  UHF,  however,  these  effects 
were  very  weak. 

C.  Spectral  Analyses 

The  spectra  of  the  complex  signal  fluctuations  were  computed  for  4-s 
periods  immediately  before  and  after  detonation.  These  spectra  were  cal¬ 
culated  using  standard  fast-Fourier-transform  techniques  after  the  (com¬ 
plex)  mean  value  of  the  data  had  been  subtracted.  Then  the  magnitudes 
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FIGURE  19  FATH  4  UHF  SIGNAL  FLUCTUATIONS 
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FIGURE  20  UHF  SIGNAL  FLUCTUATIONS  FOR  PATHS  4,  5,  AND  6 


of  the  positive  and  negative  frequency  components  were  summed,  and  the 
resulting  spectra  were  smoothed,  using  a  running  average.  Tests  with 
some  of  the  data  using  prewhitening  techniques  did  not  lead  to  any  but 
differences  in  detail,  so  these  more  sophisticated  methods  were  not  used 
for  the  most  part. 

The  fluctuation  spectra  for  the  various  frequencies  and  paths  were 
very  similar  to  one  another.  Thus  we  will  show  only  a  few  examples  here. 

Figure  21  presents  the  413-MHz  Path  1  UHF  spectra.  After  detonation,  the 

_2 

spectrum  exhibited  an  f  dependence.  The  other  Path  1  UHF  and  the  L-band 
spectra  were  virtually  identical  to  this  one,  and  the  S-band  spectrum  was 
very  similar.  At  X-band,  however,  there  was  less  difference  between  the 
pre-  and  post-detonation  spectra,  and  the  spectrum  was  flatter  below  about 
20  Hz  (Figure  22).  There  also  appears  to  be  a  peak  between  20  and  50  Hz, 
which  also  occurred  to  a  lesser  degree  at  S-band.  This  peak  also  appears 
in  the  Path  3  X-band  spectrum  (Figure  23),  which  is  virtually  identical 
to  the  Path  1  and  Path  2  spectra  as  well.  In  general,  except  for  being 
slightly  weaker,  the  post-detonation  spectra  for  the  paths  offset  from 
ground  zero  were  the  same  as  the  ones  at  the  same  wavelength  for  Path  1. 
Furthermore,  Figure  24  shows  that  the  Path  4  X-band  spectrum  is  only 
slightly  different  from  the  Path  1  spectrum  (Figure  22),  even  though  the 
dust  cloud  had  not  yet  reached  Path  4.  This  result  may  show  that  turbu¬ 
lent  air  density  fluctuations  contributed  significantly  to  signal  fluctua¬ 
tions  for  paths  penetrating  the  dust  cloud. 


IV  INTERPRETATION 

Figure  25  compares  the  measured  413-MHz  path  phase  shift  to  that  com¬ 
puted  using  a  simple  model  for  the  average  dust  density  and  the  measured 
width  of  the  cloud.  As  discussed  in  Appendix  B,  this  model  is  based  on 
the  measured  volume  of  the  DICE  THROW  dust  cloud.  Figure  26,  which  is 
the  same  as  Figure  B-38,  presents  the  average  dust  density  as  a  function 
of  time.  At  first,  it  is  assumed  that  10  kt  of  soil  was  aloft.  As  time 
went  on,  the  amount  of  soil  aloft  is  assumed  to  decrease  so  as  to  cause 
the  average  dust  density  to  follow  the  interrupted  line  in  Figure  26. 

After  T  +  20  s,  the  amount  of  mass  aloft  is  assumed  to  have  been  constant, 
and  further  average  dust  density  decreases  were  the  result  of  cloud  volume 
increases.  The  computed  absorption  is  also  compared  to  the  measured  ampli¬ 
tude  fluctuations.  The  averages  of  the  properties  of  the  soil  samples 
summarized  in  Table  A-3  were  used  for  these  computations.  Before  T  +  4  s, 
this  model  only  predicts  about  one-third  of  the  measured  phase  shift.  The 
discrepancy  can  easily  be  accounted  for  by  nonuniformity  of  dust  density 
within  the  cloud. 

At  these  early  times,  the  dust  density  was  certainly  much  greater  in 
the  center  of  the  cloud.  The  much  smaller  phase  shifts  on  Paths  2  and  3 
(Figure  12)  demonstrate  that.  Also,  the  comparison  made  in  Figure  B-4 
shows  that  the  crater  was  much  smaller  in  extent  than  the  dust  cloud. 
Greater  dust  densities  would  be  expected  nt«r  the  crater.  Consider  a 
simple,  cylindrical  model  of  the  dust  cloud  In  which  the  dust  density  de¬ 
creases  exponentially  with  radius  from  the  center  to  the  edge.  For  the 
same  average  density  as  with  the  uniform  cloud,  the  peak  density  at  the 
center  would  only  have  to  be  four  times  the  average  value  to  increase  the 
phase  shift  by  a  factor  of  three.  At  the  boundary,  the  density  for  this 
nonuniform  model  would  be  about  one-third  of  the  average  value.  This 
arbitrary  model  illustrates  how  an  actually  rather  modest  amount  of  non¬ 
uniformity  can  greatly  affect  the  result. 
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i  FIGURE  25  COMPARISONS  BETWEEN  COMPUTED  ABSORPTION  AND  PHASE 

SHIFT  AND  MEASURED  AMPLITUDE  AND  PHASE  CHANGES 
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At  413  MHz,  absorption  was  clearly  negligible  compared  to  the  effects 
of  diffractive  fading.  Nor  is  it  likely  that  scattering  losses  were  sig¬ 
nificant  either.  Of  course,  these  conclusions  apply  for  the  small  DICE 
THROW  event,  and  it  should  be  recognized  that  UHF  systems  are  not  auto¬ 
matically  immune  from  degradations  due  to  dust  clouds. 

Distortions  caused  by  the  very  complex  diffraction  phenomena  make  it 
very  difficult  to  determine  the  time-phase-shif t-versus-time  behavior. 

The  interrupted  curve  in  Figure  2 5  represents  our  best  estimate  of  the 
413-MHz  phase  shift.  We  believe  that  that  signal  underwent  a  phase  shift 
larger  than  that  measured  before  T  +  2  s  because  of  the  positive  phase 
shift  seen  on  Path  2  (Figure  12),  where  the  fading  pattern  was  very  simi¬ 
lar  co  that  of  Path  1. 

The  interrupted  curve  in  Figure  25  is  the  basis  for  computing  the 
absorption  at  higher  frequencies.  Figure  27  shows  comparisons  between 
the  computed  absorption  and  the  actual  signal  power  fluctuations  measured 
at  L-,  S-,  and  X-band.  At  all  three  frequencies,  che  computed  absorption 
was  less  than  the  measured  signal-strength  decline.  The  rest  of  the  loss 
was  caused  by  scattering.  At  L-band,  the  total  extinction  curve  appears 
to  have  roughly  the  same  shape  as  the  calculated  absorption  curve.  Even 
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FIGURE  26  AVERAGE  DUST  DENSITY  IN  DICE  THROW  CLOUD  INFERRED 
FROM  VOLUME  MEASUREMENTS 


though  the  coherent  signal  was  lost  during  part  of  the  time,  the  scatter¬ 
ing  loss,  which  was  approximately  equal  to  the  absorption  loss,  was  pro¬ 
portional  to  the  integrated  dust  density  along  the  signal  path. 

That  was  not  the  case  at  S-  and  X-band.  Both  of  those  signals  ex¬ 
perienced  a  sort  of  saturation  effect,  whereby  the  signal  amplitudes  fell 
to  low,  but  steady  (except  for  short-term  fluctuations)  values  until 
rather  rapid  recoveries  occurred.  Note  that  the  recovery  at  S-band  sig¬ 
nificantly  preceded  that  at  X-band.  This  means  that  the  recovery  was  not 
caused  by  a  simple  decrease  of  material  in  the  signal  path.  Unlike  the 
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L-band  case,  multiple  scattering  evidently  was  very  important  at  the  two 
higher  frequencies.  One  of  the  effects  of  multiple  scattering  would  be 
to  redirect  energy  back  into  the  forward  direction,  thus  limiting  the 
amount  of  scattering  loss.  Multiple  scattering  is  also  apparently  the 
explanation  why  the  measured  extinction  in  decibels  was  not  approximately 
proportional  to  frequency.  Based  on  the  S-band  result,  there  should  have 
been  a  45-dB  signal  decrease  at  X-band. 

As  the  dust  density  decreased  and  larger  particles  fell  out,  a  criti¬ 
cal  value  of  density  or  of  the  size-distribution  parameter,  Y>  was  reached 
and  the  extinction  decreased  dramatically.  Because  the  integrated  dust 
density  was  changing  very  gradually,  and,  as  discussed  in  connection  with 
Figures  A-6  and  A-7,  because  of  the  extreme  sensitivity  of  the  scattering 
loss  to  the  value  of  y,  this  threshold  was  probably  primarily  controlled 
by  the  size-distribution  parameter. 

In  Section  III-B  we  noted  that  after  initial  fallout,  the  phase 

shifts  (see  Figure  17)  corresponded  to  an  integrated  dust  density  of 

2 

2.7  g/cm  .  This  value  stayed  rather  constant  until  the  dust  cloud  was 

blown  out  of  the  signal  path.  During  this  time  interval  (i.e.,  T  +  10  s 

to  T  +  20  s),  the  cloud  was  about  560  m  across  (Figure  B-34),  giving  an 
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average  dust  density  of  4.8  x  10  g/cm  .  The  Path  4  occultation  (see 

Figure  18)  gave  rise  to  a  phase  shift  that  corresponded  to  a  maximum 

2 

integrated  dust  density  of  about  3.5  g/cm  .  However,  because  the  signal 

path  penetrated  only  a  small  part  of  the  cloud,  the  maximum  affected  path 

length  was  only  about  125  m.  Thus  the  average  dust  density  in  that  upper 
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part  of  the  cloud  was  approximately  2.8  X  10  g/cm  .  After  the  initial 
fallout  phase,  the  dust  density  apparently  was  much  greater  on  the  upper 
parts  of  the  cloud  than  toward  the  base.  Our  crude  model  based  on  mea¬ 
sured  cloud  volume  and  estimates  of  soil  initially  lofted  and  of  that 
eventually  missing  from  the  crater  predicted  dust  densities  (Figures  26 
and  B-38)  that  are  in  reasonably  good  agreement  with  these  values. 
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A  dust  cloud  density  in  the  range  of  10  to  10  g/cm  is  also 
consistent  with  observed  attenuation  and  computed  extinction.  If 
multiple-scattering  phenomena  are  ignored  and  all  of  the  attenuation  is 
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accounted  for  by  scattering  loss,  the  extinction  coefficients  at  S-  aid 

X-bands  are  roughly  7  and  10  dB/km  respectively.  According  to  the  results 

shown  in  Figure  A-6  for  our  simple  extinction  model,  the  average  cloud 
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mass  density  would  only  have  to  be  about  10  g/cm  for  y  =  3.5.  Because 
of  the  effects  of  multiple  scattering,  this  can  be  considered  an  upper 
limit,  within  the  accuracy  of  the  model.  This  means  that  large  particles 
comprising  only  a  small  fraction  of  the  mass  in  the  dust  cloud  can  easily 
account  for  the  measured  extinction,  giving  a  consistent  interpretation 
of  the  observed  behavior. 

As  the  larger  particles  fell  out,  the  effective  value  of  V  decreased 
until  scattering  losses  suddenly  diminished.  However,  the  total  mass  in 
the  dust  cloud  apparently  declined  only  a  little  during  that  time.  Accord¬ 
ing  to  Figure  A-5,  the  X-band  result  shows  that  there  were  a  substantial 
number  of  particles  larger  than  0.1  cm  in  the  signal  path  until  about 
T  +  3  s.  Since  the  critical  particle  size  at  L-band  was  about  1  cm,  there 
were  numerous  particles  larger  than  that  present. 

In  general,  except  when  the  threshold  combination  of  density  and  Y 
was  exceeded,  rf  effects  were  fairly  mild.  Most  of  that  can  be  ascribed 
to  the  small  size  of  the  DICE  THROW  event.  In  the  nuclear  situation,  how¬ 
ever,  where  much  larger  amounts  of  soil  can  be  lofted  over  much  larger 
areas  and  volumes  for  much  longer  periods  of  time,  the  cumulative  effect 
of  even  low-density  dust  clouds  over  long  signal  paths  could  be  very 
serious.  And,  once  the  threshold  is  crossed,  a  signal  can  deteriorate 
very  quickly.  It  is  probable  that  the  fallout  of  very  likely  less  than 
a  few  percent  of  the  total  lofted  mass  during  DICE  THROW  effectively  led 
to  signal  recovery  at  the  higher  frequencies.  Conversely,  a  patch  of 
particles  within  a  much  larger  cloud  having  characteristics  exceeding  the 
threshold  could  lead  to  a  decline  in  signal  strength  and  quality  even 
though  the  density  changed  very  little.  Furthermore,  the  loss  of  a  coherent 
signal  component  would  cause  a  very  serious  decline  in  signal  quality  for 
many  systems  that  could  otherwise  compensate  for  large  changes  in  signal 
levels. 
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V  CONCLUSION 


The  DICE  THROW  UHF/SHF  coherent-transmission  experiment  was  generally 
very  successful.  Although  three  signals  were  lost  after  the  detonation 
due  to  shockwave  effects,  all  14  signal  channels  functioned  properly  for 
at  least  the  first  10  seconds,  during  which  the  most  significant  rf 
effects  occurred.  There  is  fairly  good  agreement  between  the  computed 
effects  using  a  very  simple  model  of  the  dust  cloud  and  the  measured  sig¬ 
nal  perturbations.  A  consistent  explanation  of  the  effects  is  thereby 
obtained. 

For  much  of  the  time,  the  DICE  THROW  rf  effects  were  very  mild. 

Most  of  this  can  be  ascribed  to  the  small  size  of  the  detonation.  Besides 
that,  the  capped-cylinder  ANFO  charge  configuration  also  tends  to  produce 
low  flow-field  velocities,  resulting  in  rapid  fallout  and  early  cessation 
of  strong  effects.  Worse  yet,  the  downdraft  in  the  center  of  the  flow 
field  intensifies  the  fallout  rate  of  large  particles.  Even  at  the  same 
yield,  in  the  nuclear  situation  the  higher  velocities  and  the  updraft 
would  keep  larger  particles  and  more  overall  mass  aloft  for  much  longer 
periods  of  time.  Strong  rf  degradations  and  perturbations  would,  there¬ 
fore,  be  expected  to  persist  much  longer  than  simple  yield-scaling  would 
predict.  Scaling  from  the  DICE  THROW  configuration  to  the  nuclear  situa¬ 
tion  needs  to  be  quantified. 

The  onset  of  and,  more  significantly,  the  decay  of  strong  effects 
were  very  sudden.  There  appears  to  have  been  a  threshold  combination  of 
dust  density  and  particle-size  distribution.  Once  the  threshold  was 
crossed,  average  extinction  coefficients  rose  to  as  much  as  10  dB/km. 
Locally,  this  attenuation  must  have  been  much  larger,  perhaps  reaching 
30  dB/km  or  more.  This  was  accompanied  by  complete  losses  of  the  coherent 
signal  components,  and  apparently,  multiple  scattering  became  dominant. 
Besides  the  large  decrease  in  average  signal  strength,  very  strong  fading 


also  occurred,  which  was  more  variable  than  Rayleigh  fading.  All  of 
these  effects  were  frequency  dependent,  becoming  worse  as  the  frequency 
increased. 
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Inferred  average  dust  densities  ranged  between  10  and  10  g/cm  , 
depending  on  the  part  of  the  DICE  THROW  cloud  in  the  signal  path  and  on 
the  time  after  detonation.  However,  even  though  the  average  density 
along  a  signal  path  was  low,  nonuniformity  apparently  caused  the  threshold 
condition  to  be  crossed  locally,  resulting  in  very  large  losses  from  the 
higher-density  patches.  It  would  be  likely  that  in  the  nuclear  case  the 
threshold  would  be  exceeded  for  much  longer  periods  of  time  and  that 
scattering  losses  would  be  similar  to  those  responsible  for  the  extinction 
coefficients  measured  during  DICE  THROW,  The  possibility  of  encountering 
local  high-loss  patches  within  a  cloud  should  be  considered  when  evaluat¬ 
ing  the  performance  of  any  given  system.  Attenuation  in  such  patches 
could  be  as  great  as  30  dB/km  at  S-band  and  above. 

Although  they  are  not  as  dramatic,  the  effects  of  a  low-density  or 

relatively  homogeneous  cloud  of  sufficient  extent  must  be  considered  as 

well.  In  that  case,  the  attenuation,  which  is  dominated  by  absorption, 

could  become  sufficiently  large  to  affect  certain  systems.  Even  at  UHF, 

where  the  effects  of  dust  clouds  are  fairly  weak,  the  attenuation  could 

be  greater  than  tolerable  for  a  low-margin  system.  For  example,  suppose 

there  were  a  30-km-wide  region  of  dust  having  an  average  density  of 
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3  x  10  g/cm  .  Then,  using  the  average  dielectric  properties  of  the 
DICE  THROW  soil  samples,  the  absorption  would  rise  from  a  possibly  sig¬ 
nificant  1.4  dB  at  300  MHz  to  a  definitely  serious  value  of  14  dB  at 
10  GHz. 

Unlike  scattering  losses,  absorption  is  proportional  to  the  loss 
tangent  of  the  material  comprising  the  dust  particles,  and  is  therefore 
very  sensitive  to  moisture  content.  In  the  nuclear  situation,  then,  the 
effects  of  heating  by  the  hot  fireball  on  the  soil  lofted  needs  to  be 
considered.  The  amounts  of  water  vaporized  and  recondensed  are  the  im¬ 
portant  questions.  If  a  large  amount  of  water  initially  locked  up  in 
crystalline  structures  were  liberated,  it  is  conceivable  that  absorption 


losses  could  be' greater  after  recondensation  than  for  the  same  cloud 
composed  of  unheated  particles.  Although  predicted  temperatures  were 
several  hundred  degrees  Celsius  above  ambient  in  the  DICE  THROW  cloud  for 
a  number  of  seconds,  heating  effects  were  undoubtedly  not  as  significant 
as  they  would  have  been  for  a  nuclear  detonation  of  the  same  yield. 
Nevertheless,  heating  of  clear  air  outside  the  dust-laden  region  in  the 
DICE  THRCW  flow  field  resulted  in  a  phase  shift  comparable  in  magnitude 
to  that  produced  by  the  dust.  This  effect  tended  to  enhance  refractivity 
gradients,  which  in  turn  would  lead  to  more  fading  and  dispersion. 

A  test  such  as  DICE  THRCW  only  reproduces  some  aspects  of  the  nuclear 
environment.  In  general,  rf  effects  were  weaker  than  would  be  expected 
from  a  nuclear  event  that  lofted  the  same  amount  of  soil  initially.  How¬ 
ever,  the  results  from  the  DICE  THROW  UHF/SHF  coherent  transmission  experi¬ 
ment  showed  that  very  serious  signal  degradations  can  occur  and  that  there 
are  a  number  of  interacting  phenomena  that  need  to  be  considered. 


Appendix  A 


EFFECTS  OF  DUST  CLOUDS 
ON  ELECTROMAGNETIC  PROPAGATION 

1.  Introduction 

Figure  A-l  summarizes  the  phenomena  associated  with  the  interaction 
of  electromagnetic  waves  with  dust  clouds.  Part  of  the  incident  energy 
is  absorbed  by  the  dust  cloud.  Matter  in  the  dust  cloud  also  causes  dif- 
ferential  phase  retardations,  causing  irregular  emergent  wave  fronts, 
which  in  turn  lead  to  angle~of -arrival  fluctuations  and  diffractive  phe¬ 
nomena  associated  with  propagation  over  distance  that  manifest  themselves 
as  scintillations  or  fading  behavior.  Frequency-selective  fading  (coherent- 
channel-bandwidth  reduction)  can  also  occur.  Under  most  circumstances 
(e.g.,  when  all  the  dust  particles  are  small  compared  to  the  wavelength), 
absorption  and  phase  shifts  are  thought  to  be  closely  related  aspects  of 
the  same  phenomenon--a  change  in  the  bulk  complex  dielectric  properties 
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FIGURE  A-1  RF  PROPAGATION  AND  DUST  CLOUDS 
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of  the  propagation  medium  (dust-laden  air)  that  can  be  related  to  the 
mass-density  and  electromagnetic  properties  of  the  individual  solid  grains. 
Under  the  small-particle  approximation,  the  particle  size  distribution  does 
not  affect  the  results. 

Besides  the  absorption  loss,  scattering  losses  also  contribute  to  the 
extinction,  or  attenuation,  of  the  electromagnetic  waves.  Unlike  absorp¬ 
tion,  which  is  usually  treated  as  a  coherent,  bulk  (or  "fluid")  phenomenon, 
scattering  is  considered  an  incoherent  or  semicoherent  phenomenon  that 
depends  crucially  either  on  the  particle  size  distribution,  the  spectrum 
of  (turbulent)  dust  density  fluctuations,  or  both. 

Backscatter,  which  can  cause  radar  clutter  interference,  is  only  one 
aspect  of  the  scattering  problem.  Indeed,  all  the  phenomena  mentioned 
above  are  merely  different  aspects  of  the  same  thing — the  interaction 
of  the  electromagnetic  energy  with  matter  in  a  grain  of  dust.  However, 
rather  than  trying  to  treat  this  interaction  in  a  general  way,  it  has 
been  far  more  convenient  to  study  its  different  aspects  in  more  restric¬ 
tive  ways  that  try  to  isolate  the  most  important  and  relevant  parameters 
for  the  application.  The  division  into  coherent  and  incoherent  or  semico¬ 
herent  effects  mentioned  above  is  an  example. 

Several  aspects  of  the  interaction  of  electromagnetic  waves  with  dust 
clouds  have  been  reported  on  by  SRI  International  as  part  of  this  or 
closely  related  work.  Bollen  et  al.6  considered  the  bulk  phase  shift  and 
scattering  loss  problem  in  the  DICE  THROW  effects-predictions  part  of 
their  discussion  of  dust  clouds  and  crater  ejecta  from  nuclear  and  HE 
detonations.  Burns  and  Winkleman3  related  measured  solid-grain  dielectric 
and  absorptive  properties  to  those  of  the  DICE  THROW  Main  Event  dust  cloud. 
Chesnut6  investigated  the  relative  intensities  of  incoherent  backscattering 
and  backscattering  from  density  fluctuations  in  turbulent  dust  clouds. 

This  appendix  is  primarily  a  summary  of  the  results  in  those  reports, 
which  will  be  used  to  interpret  the  experimental  results. 


Bulk  Absorption  and  Phase  Shift 


The  relationship  between  the  dielectric  properties  of  powders  (a  dust 
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cloud  is  a  very-low-density  powder)  and  those  of  the  solid  grains  is 
specified  by  a  so-called  "mixing  formula,"  Some  experimental  evidence7 
supports  the  use  of  the  Rayleigh  mixing  formula, 


±  It ili\ .  ± 

M‘d  +  2/  M‘.  +  2/ 


(A-l) 


where 


P^  *  Mass  density  of  powdered  material 

P  *  Mass  density  of  solid  grains 
s 

*  Complex  relative  dielectric  constant  of  powder 

eg  *  Complex  relative  dielectric  constant  of  solid. 

This  formula,  which  is  closely  connected  to  the  Lorentz-Lorenz  or  Clausius- 
Massotti  equations,8’9  is  based  on  the  approximation  that  all  the  grains 
have  induced  dipole  moments.  This  means  that  they  are  all  smaller  than 
the  wavelength  and  are  Rayleigh  scatterers.  The  forward  (i.e. ,  in  the 
direction  of  propagation)  scattered  waves  add  up  coherently  according  to 
the  Ewald-Oseen  extinction  theorem  in  such  a  way  that  the  properties  of 
the  resulting  wave  can  be  obtained  using  Eq.  (A-l). 

The  resulting  absorption,  a,  and  phase  constants,  g,  for  waves  in 


dust  clouds  are:v 


2  \pJ(sr  +  2)2 

e“4+2(3M 


where 


(A-2a) 


(A- 2b) 


k  =  2tt/\  =  Wave  number 


e  =  Real  part  of  e 
r  r  s 

tan  6  =  e$/er  =  Loss  tangent 
=  Imaginary  part  of  eg 
p^/pg  =  Relative  dust  cloud  density. 


The  intrinsic  loss  terra,  tan  6q,  and  the  conductivity,  Oq ,  and  e^, 

were  found  to  be  nearly  independent  of  frequency  from  0,1  to  12  GHz 
9  -1 

[eQ  »  (36tt  X  10  )  in  MKS  units].  Table  A-l  summarizes  the  measurements 
made  on  the  three  DICE  THROW  soil  samples  collected  after  the  event  from 
the  crater  and  lip  ("sand,"  "caliche  A,"  and  "caliche  B").  While  did 
not  change  significantly  as  the  samples  were  dried,  both  tan  6q  and  oQ 
were  very  sensitive  functions  of  moisture  content.  The  solid-grain  prop¬ 
erties  inferred  using  the  Rayleigh  mixing  law  are  given  in  Table  A-2  for 
1  and  10  GHz. 


Equation  (A-5) ,  and  some  equations  leading  to  the  equivalent  of  Eq.  (A-2a)  , 
were  in  error  in  Ref.  3;  the  conclusions  drawn  there  were  unaffected,  be¬ 
cause  the  effects  of  the  errors  were  slight. 
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Table  A- l 


SUMMARY  OF  DIELECTRIC  PROPERTY  MEASUREMENTS 


Sample 

Moisture 

Content 

(percent) 

Density 

in 

«r 

tan  6 

0 

°o 

Fixture 

(g/cm3) 

1  GHz 

10  GHz 

(mmho/m) 

Sand  (undried) 

2.74 

1.521 

2.9 

2.5 

0.025 

2.4 

Sand  (dried) t 

1.512 

2.7 

2.5 

0.004 

<5  xl0“3 

Sand  (dried  and 
exposed) * 

0.90 

1.549 

2.8 

2.5 

0.008 

<10-2 

Caliche  A  (undried) 

12.73 

1.207 

2.8 

2.4 

0.024 

3.1 

Caliche  A  (dried 
and  exposed) 

3.38 

1.200 

2.5 

2.1 

0.00? 

<10-2 

Caliche  A  (dried) 

w0  0 

1.201 

2.4 

2.1 

0.006 

<10"2 

Caliche  B  (undried) 

9.03 

1.472 

3.0 

2.7 

0.021 

1.4 

Caliche  B  (dried) 

1.388 

2.6 

2.2 

0.007 

Cs! 

• 

O 

*~-4 

V 

As  collected. 

^ Baked  for  A  days  in  vacuum  chamber  at  U0°C. 

* 

Exposed  to  ambient  conditions  for  several  days  after  being  dried. 


Table  A-2 

INFERRED  AVERAGE  DIELECTRIC  CONSTANT,  MOLECULAR  LOSS  TANGENT, 
AND  CONDUCTIVITY  FOR  SOLID  PARTICLES 


Solid-Particle 
Density 
(g/ cm-*) 

Moisture 

Content 

(percent) 

f 

(GHz) 

Solid-Particle 

Properties 

€ 

r 

tan  6 

0 

°o 

(mmho/m) 

Sand 

2.56 

2.74 

mm 

6.6 

0.051 

m 

99 

4.8 

0.043 

mm 

Caliche  A 

2,64 

12.73 

i 

14.7 

H 

10 

7.9 

■Sfl 

o3 


i 

t 

; 


£ 


The  curves  shown  in  Figure  A-2  result  from  inserting  the  parameter 

values  listed  in  Table  A-2  into  Eq.  (A-5).  The  conductivity  term  does 

not  appreciably  contribute  to  absorption  above  about  1  GHz.  Table  A-3 

lists  the  excess  phase  retardation  and  absorption  per  kilometer  normalized 
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to  a  uniform  dust  cloud  density  of  10  g/cm  at  several  frequencies  in 
the  0.3-to-10-GHz  range  for  the  DICE  THROW  samples.  Finally,  Figure  A-3 
presents  the  excess  phase  retardation  and  absorption  per  kilometer  as  a 
function  of  uniform  dust  density  computed  at  1  and  10  GHz.  The  dielectric 
property  measurements  and  calculations  presented  in  this  appendix  were 
used  to  interpret  the  experimental  results. 

3.  Scattering  Loss 

The  other  component  of  extinction  is  loss  due  to  energy  scattered 
away  from  the  forward  direction.  Two  processes  are  responsible  for  the 
scattering.  The  first  process  is  incoherent  scattering,  where  the  power 
scattered  from  each  grain  is  summed  for  all  of  the  grains.  Chesnut6  has 
termed  the  other  process  "semicoherent."  This  type  of  scattering  arises 
from  nonuniformities  in  the  particulate  cloud  density,  and  depends  on  the 
strength  and  scale  size  of  the  density  fluctuations.  It  is  termed  "semico¬ 
herent"  because  scattered  fields  from  nearby  particles  are  added  coherently. 
A  limiting  case  is  the  example  of  a  perfectly  uniform  dust  cloud  where  the 
only  density  discontinuity  occurs  at  its  edge;  then  the  scattering  can  be 
found  using  the  results  of  the  preceding  section  to  determine  boundary 
reflection  coefficients.  In  that  example,  scattering  from  the  grains  acts 
in  a  completely  coherent  way  within  the  cloud.  However,  very  abrupt  bound¬ 
aries  (i.e.  ,  transition  distances  small  compared  to  the  rf  wavelength)  are 
required  to  cause  appreciable  scattering.  In  general,  incoherent  scatter¬ 
ing  has  been  consistently  considered  the  dominant  process  of  scattering 
from  particulate  clouds. 

Because  incoherent-scatter  strength  depends  strongly  on  the  particle- 
size-aistribucion  function,  it  is  necessary  to  adopt  a  model  for  the  grain 
size  function  in  a  dust  cloud.  Although  there  is  very  little  or  no  theore¬ 
tical  justification  for  such  a  choice,  an  inverse-power-law  function  is 
apparently  often  picked  on  the  basis  of  empirical,  soil-sieving  evidence. 
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FIGURE  A-3  EXCESS  PHASE  RETARDATION  AND  ABSORPTION  PER 
KILOMETER 


ABSORPTION  —  dB/km 


As  an  example,  this  is  the  choice  made  by  Kownacki  and  Welty,10  and  it  is 
the  choice  made  here,  although  we  will  modify  it  somewhat.  Then,  if  a  is 
the  grain  size  (e.g. ,  the  radius,  for  spherical  grains), 


n(a)  =  Ka 


(A- 7) 


is  the  number  of  particles  per  unit  volume  between  a  and  a  +  da.  Typical 
values  for  y  range  from  3.5  to  9. 11  For  spherical  grains  the  mass  density 
in  the  dust  cloud  is  given  by 
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(A-8) 


For  a  particle  in  the  Rayleigh  regime  (where  ka  «  1) ,  the  scattering 
and  absorption  cross  sections  are  given,  respectively,  by 
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(A- 10) 


when  tan  6  «  1. 

When  ka  «  1,  a  »  a  and  extinction  is  essentially  all  due  to  ab- 
3  S 

sorption.  The  ratio  of  wavelength  to  radius  for  which  o  =  a  is 
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X  2(er  •  X) 

a  =  2n  9  er  tan  6 


(A-ll) 


For  our  typical  soil  values — say,  ■  5  and  tan  6  =  0.03 — X/a  =  18.  This 
value  is  only  about  one-third  that  where  the  Rayleigh  and  optical-scattering 
asymptotes  intersect  (Figure  A-4) .  We  therefore  conclude  that  scattering 
loss  is  negligible  for  the  relevant  materials  here  ir  the  range  where  the 
Rayleigh  approximation  is  valid. 

Absorption  by  grains  in  the  Rayleigh  region  was  considered  in  the 
preceding  section.  Because 
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FIGURE  A-4  INTERSECTIONS  OF  OPTICAL  AND  RAYLEIGH  ABSORPTION  AND  SCATTERING 
ASYMPTOTES 


69 


a  *  1/2  a(^  =  1/2  /  a„n(a)da  - 


/*  6nk  e  tan6  r 

a  n(a)da  - - - - s-  /  i 

U  +  2)1  J 


a3n(a)da  (A-12)* 


and,  from  Eq.  (A-8) , 
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Eq.  (A-2a)  follows.  Not  surprisingly,  the  same  expression  for  absorption 
can  be  derived  using  either  approach. 


Extinction  when  all  of  the  grains  are  Rayleigh  scatterers  seems, 
therefore,  to  be  well  understood  (semicoherent  scattering--to  be  considered 
below--is  probably  much  weaker  than  the  incoherent  process).4  On  the 
other  hand,  the  situation  is  very  uncertain  when  many  particles  are  in 
the  optical  region.  Mie  resonances  may  occur,  may  not  occur,  or  may 
be  significant,  depending  on  the  geometric  regularities  (e.g. ,  perfect 
sphericity)  of  the  particles  and  the  loss  tangent  of  the  material  compris¬ 
ing  the  scatterers.  Also,  there  has  been  some  question5*6  whether  the 
extinction  cross  section  in  the  optical  regime  is  twice  the  geometrical 
cross  section  oi  equal  to  the  geometrical  cross  section. 

Q 

Born  and  Wolf  assert  that  the  extinction  cross  section  is  twice  the 
geometrical  cross  section  for  large,  totally  absorbing  obstacles.  They 
also  report  the  work  of  other  authors  that  support  this  conclusion  for 
perfectly  reflecting,  nonabsorbing,  and  weakly  absorbing  spheres.  Other 
calculations10  provide  further  support  for  that  result.  The  question 
seems  to  be  whether  or  not  energy  diffracted  in  the  forward  direction  into 
the  region  around  the  geometric  shadow  is  truly  "lost"  from  the  incident 
wave  and  should  be  counted  as  extinction.  In  this  report,  we  will  assume 
that  the  short-wavelength  extinction  cross  section  oQ  is  twice  the  geome¬ 
tric  one,  and,  for  simplicity,  essentially  ignore  any  possible  Mie  reso¬ 
nances.  Our  extinction  model  will  then  be  the  one  shown  in  Figure  A-4 


The  limits  of  integration  in  Eq.  (A-12)  should  be  over  all  sizes,  if  all 
of  the  grains  are  Rayleigh  scatterers,  or,  if  there  are  large  particles, 
up  to  a  limit  (ka  —  0.25),  obtained  from  considerations  discussed  below. 
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by  the  dashed  line,  which  is  twice  the  scattering  cross  section,  for 
ka  £  0.25  (o  was  set  to  2  a  on  the  Rayleigh  asymptcte  to  compensate 

6  S 

roughly  for  any  Mie  resonances,  which  tend  to  raise  og  in  the  0.25  j£  ka  s:  1 
range).  Figure  A-5  shows  the  values  of  a  for  which  ka  =  0,25  and  ka  *  1, 
which  turns  out  to  be  very  close  to  the  value  (for  er  ■  5)  where  the  Ray¬ 
leigh  and  optical  asymptotes  intersect,  for  the  frequency  interval  rele¬ 
vant  to  the  DICE  THROW  transmission  experiment.  It  is  clear  that  scatter¬ 
ing  losses  could  have  been  very  important  for  the  higher  measurement 
frequencies  at  early  times  after  detonation,  when  many  centimeter  and 
larger  particles  were  probably  airborne. 

Bollen  et  al.6  made  predictions  of  DICE  THROW  transmission  experiment 
effects  using  a  particle  size  distribution  thought  to  have  been  representa¬ 
tive  of  the  DICE  THROW  site.  Their  particle  size  distribution  ("distribu¬ 
tion  A")  can  be  reasonably  well  approximated  by  a  two-segment  modification 
of  Eq.  (A-7): 
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About  half  the  mass  is  in  particles  bigger  than  5  cm.  Because  all  of 
the  particles  in  the  y  =  5  segment  (i.e. ,  those  with  a  <  0. 1  cm)  are 
below  the  ka  =  0.25  point  for  frequencies  of  interest  here  (cf.,  Figure 
A-5),  they  will  make  a  negligible  contribution  to  the  scattering  component 
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As  expected,  the  lower  integration  limit  makes  less  than  a  1%  contribution 

to  this  expression  for  or;  this  result  further  supports  ignoring  the  Y  ■  5 

segment.  The  second  term  in  Eq.  (A- 16)  comes  from  the  upper  integration 

-1/2 

limit,  and,  since  it  varies  as  a^  ,  it  becomes  relatively  less  important 
as  there  are  more  large  particles  or  as  the  frequency  increases.  Values 
for  a  obtained  using  Eq.  (A-16)  are  comparable  to  those  computed  by  Bollen 
et  al.,5  where  numerical  integration  was  used  with  the  actual  particle 
size  distribution  rather  than  with  this  model. 

Extremely  large  values  of  attenuation  can  occur  when  there  are  many 

large  scatterers  as  compared  to  absorption  by  a  cloud  of  Rayleigh  regime 

grains.  The  solid  lines  of  Figure  A-6  present  the  attenuation  computed 

3 

from  Eq.  (A-16)  at  p  =2.5  g/cm  versus  p  ,  A  comparison  with  the 

s  s 

absorption  curves  of  Figure  A-4  reveals  that  extinction  losses  for  this 
model  are  enormously  larger  than  absorption  losses  from  the  smaller  grains. 
However,  extinction  losses  are  very  sensitive  to  V.  The  dashed  line  in 
Figure  A-6  is  the  extinction  computed  for  a  model  where  y  =  5  from 


a  ■  7  x  10  m  to  a.  *  0.5  m:  this  curve  is  substantially  below  the 
s  b 


absorption  curves.  Figure  A-7  shows  that  the  ratio  of  scattering  loss 
for  V  »  5  to  that  of  Y  ■  3.5  for  these  models  is  in  the  range  of  10~^  to 


>-2 


10  over  the  frequency  interval  of  interest  here.  Larger  values  of  Y 
lead  to  even  smaller  values  for  scattering  loss. 


Chesnut  considered  the  relative  strengths  of  incoherent  aad 


semicoherent  backscattering  and  concluded  that  coherent  scattering 

was  substantially  the  weaker  for  dust  clouds.  Figures  A-8  and  A-9  sum¬ 
marize  his  results.  The  dashed  curves  apply  for  particulate  clouds  having 
only  one  grain  size,  and  the  solid  lines  apply  for  a  continuous  size  dis¬ 


tribution  dominated  by  small  particles  (Y  —  6  for  0.6  <  a  <  3  mm  and 


Y  —  9  for  a  >  3  tran).  These  results  were  also  based  on  p  ■  2g/cm  ,  an 
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outer  scale  of  100  m,  and  a  turbulent  spectrum  approximating  the  Kolmo¬ 
gorov  spectrum.  Perhaps  the  most  significant  assumption  made  was  that 
the  variance  of  particle  number  was  equal  to  its  expectation,  which 
probably  greatly  overestimates  0  0j,/^nc0h  most  circumstances.  Although 
semicoherent  scattering  is  relatively  more  important  when  particles  are 
small  and  at  low  frequencies,  it  is  very  small  in  an  absolute  sense.  Be¬ 
cause  a  preponderance  of  small  grains  is  likely  to  occur  only  at  low  dust 
cloud  mass  densities,  semicoherent  scattering  is  probably  negligible. 


4.  Summary 

Scattering  losses  are  probably  negligible  except  at  early  times 
before  large  particles  have  fallen  out  of  the  dust  cloud.  Absorption 
losses,  which  dominate  in  the  Rayleigh  regime,  can  be  estimated  from 
knowledge  of  the  dielectric  properties  of  the  solid  grains  without  requir¬ 
ing  knowledge  of  the  particle  size  distribution.  When  scattering  losses 
are  important,  however,  knowledge  of  the  particle  size  is  crucial.  Semi¬ 
coherent  scattering  due  to  density  fluctuations  is  probably  negligible 
compared  to  incoherent  scattering  for  conditions  likely  to  be  present  in 
dust  clouds.  Because  extinction  coefficients  appear  as  exponents,  it 
should  be  kept  in  mind  that  small  errors  in  estimating  those  coefficients 
can  lead  to  many  decibels  of  error  when  losses  are  appreciable. 
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FIGURE  A-8  VOLUME  SCATTERING  COEFFICIENT  vs  RADAR  FREQUENCY 
FOR  VARIOUS  OUST  PARTICLE  SIZE  DISTRIBUTIONS— 
INCOHERENT  AND  SEMICOHERENT  SCATTERING 
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FIGURE  A-9  RATIO  OF  COHERENT  TO  INCOHERENT  SCATTERING  vs.  RADAR  FREQUENCY 
FOR  VARIOUS  MODELS  OF  DUST  PARTICLE  SIZE  DISTRIBUTION 
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Appendix  B 

DICE  THROW  CLOUD  DEVELOPMENT 


1.  Introduction 

In  support  of  the  primary  rf  experiment,  photographic  coverage  of 
the  DICE  THROW  event  was  used  to  document  cloud  morphology  and  position 
versus  time  for  correlation  with  microwave  measurement  results.  Seven 
35-mm  Automax  framing  cameras  containing  color  positive  film  (Kodak  type 
5256  Ektachrome)  were  mounted  at  various  locations  along  and  roughly 
perpendicular  to  the  signal  paths.  The  camera  locations  with  respect  to 
ground  zero  (GZ)  and  the  transmission  experiment  layout  were  shown  in 
Figure  2. 

Relatively  slow  framing  rates  of  1,  2,  and  5  frames  per  second  (fps) 
were  selected.  Two  cameras  were  switched  from  5  fps  to  1  fps  after 
2  minutes  in  order  to  extend  late-time  coverage  of  cloud  development. 

Focal  lengths  of  lenses  were  varied  in  an  effort  to  increase  the  quality 
of  photographic  documentation  and  allow  for  two  different  field-oi  view 
coverages  at  each  site  (except  at  the  main  transmitter).  Table  B-l  lists 
the  camera  positions,  parameters,  and  approximate  event-time  coverages. 
Cameras  at  the  unmanned  main  transmitter  (MT)  and  phase  repeater  (PR) 
sites  were  remotely  activated  at  T  -  1  min,  whereas  the  others  were 
manually  controlled.  All  cameras  were  in  a  fixed  position  with  line  of 
sight  to  GZ,  thus  facilitating  cloud  motion  measurement. 

Photographic  documentation  of  the  DICE  THROW  event  has  provided  a 
pictorial  history  of  the  event  and  scaled  measurements  of  cloud  morphology 
in  direct  support  of  data  interpretation.  DICE  THROW  morphology  can  be 
compared  with  actual  Pre-DICE  THROW  II-2  (PDT  II-2)  morphology,  as  well 
as  to  the  predicted  ANFO  behavior  celculations  of  Pre-DICE  THROW  for 
early-,  intermediate-,  and  late-time  phenomenology,  rise  rates,  and 
velocity  fields.  Comparisons  with  PDT  II-2  are  pertinent  because  PDT  II-2 
was  used  to  predict  the  DICE  THROW  behavior,6  Other  spatial  data — i.e., 
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Main  transmitter 
Phase  repeater 
North  Oscuro  Peak 


height,  width,  and  volume  versus  time— as  well  as  a  general  description 
of  cloud  development,  aid  the  understanding  of  dust  cloud  evolution  and 
eventual  disposition. 

2.  Photographic  Time  Histories  ' 

Familiarity  with  the  overall  DICE  THROW  cloud  morphology  as  a  func¬ 
tion  of  time  is  necessary  in  order  to  correlate  microwave  transmission 
data  with  specific  cloud  features  at  the  points  of  raypath  intersections. 
In  addition,  the  macrophenomenology  of  the  DICE  THROW  event,  a  628- ton 
ANFO  detonation,  can  be  a  useful  form  of  comparison  between  HE  and  nuclear 
detonations,  and  among  different  HE  charge  configurations.  The  following 
description  of  explosive  geometries  and  cloud  motion  of  DICE  THROW  is 
divided  into  time  domains  designated  as  early,  intermediate,  and  late. 

The  time  coverages  selected  are  deemed  appropriate  for  an  HE  detonation 
produced  by  an  explosive  weight  of  hundreds  of  tons,6  Composites  of 
scaled  photographs  exposed  by  cameras  at  Sites  R-i,  PR,  NOP,  and  MT  (see 
Figure  2)  are  presented  together  according  to  times  after  detonation. 

Where  appropriate,  photographs  of  HE  and  nuclear  detonations  are  presented 
to  support  comparative  observations. 

3.  Early  Time— T  to  T  +  6  s 

o  _ 

Figure  B-l  is  a  pro-detonation  view  of  the  charge  at  GZ  looking  from 
the  phase  repeater  site.  The  ANFO  charge  was  11, A  m  high  from  ground 
level  to  the  top  of  the  cylinder  cap  and  9. 1  m  in  diameter.  Figures  B-2 
through  B-6  show  early-time  (0  to  6  seconds)  view  of  DICE  THROW  from  as 
many  as  four  camera  sites  (R-l,  PR,  NOP,  and  WT).  In  Figure  B-2(a),  at 
T  ♦  0. 1  s,  heated  detonation  products  are  evident  through  the  lower  two- 
thirds  of  the  cloud,  This  material  had  the  typical  white  appearance  of 
an  ANFO  detonation  at  first  (Figures  B-2 (hi  and  (c),  T  +  0.2  si  and 
became  dome-shaped  as  it  rose.  This  mound  of  material  was  most  likely  a 
consequence  of  shockwave  reflection  at  the  surface;  that  is,  the  upward 
velocity  of  material  is  greatest  directly  over  the  detonation  point  since 
initial  upward  velocity  is  proportional  to  the  pressure  of  the  direct 


FIGURE  B-1  DICE  THROW  CAPPED  CYLINDER  CHARGE  OF  ANFO— PRE-ZERO  VIEW  FROM  PHASE 
REPEATER  SITE 


FIGURE  B-2  DICE  THROW  EARLY  TIMES  T  +  0.1  $  TO  T  +  0.5  s,  FROM  FOUR  CAMERA  SITES 


FIGURE  B-3  OICE  THROW  WILSON  CLOUD  AT  T  +  1  *.  VIEWED  FROM  PHASE  REPEATER 


shockwave.  By  T  +  1  s,  the  cloud  darkened  in  color,  presumably  due  to  the 
presence  of  unburned  carbon  from  both  the  sacking  material  and  the  non- 
stoichiometric  charge  composition  (there  was  a  slight  overabundance  of 
fuel  oil).  In  contrast  to  DICE  THROW,  no  dark  material  was  observed 
during  PDT  II-2,  apparently  due  to  obscuration  by  the  spires  of  heavy 
ejecta  and  base  surge  (Figure  B-7(a)l;  however,  such  material  was  evident 
in  the  Canadian  hemispherical  100- tort  ANF0  detonation  in  1969  (Figure  B-7 
(b)].3  A  small  Wilson  (condensation)  cloud  lasted  from  T  +  0.5  s  to 
T  +  1.2  s  (Figure  B-3). 

An  overhead  view  (Figure  B-8)  at  T  +  0.8  s  (with  the  apparent  crater 
rim  outline  approximately  drawn  to  scale)  reveals  some  interesting  features 
Dust  billows  radiating  out  from  the  cloud  center  like  wheel  spokes  are 
aligned  along  existing  dirt  roads  near  GZ.  Vehicle  traffic  churned  the 
dirt  along  these  roads  and  the  dirt  was  looser  to  greater  depths  than  the 
soil  on  surrounding  surfaces.  This  dust  was  more  readily  lofted  than 
that  from  other,  less  disturbed  areas. 

The  early-time  morphology  of  DICE  THROW  also  differed  from  that  of 
PDT  II-2  in  another  respect.  In  PDT  II-2,  large,  water-saturated  clay 
"missiles"  were  radially  ejected  or  thrown  out  from  the  center  at  high 
angles,  forming  very  noticeable  spires.  This  effect  was  not  characteris¬ 
tic  of  DICE  THROW,  apparently  because  the  water  table  at  the  DICE  THROW 
Main  Event  site  was  at  least  25  m  deep,  compared  to  only  2-3  m  at  the 
PDT  site.  However,  the  conditions  near  the  surface  appeared  similar  for 
both  events.  Therefore,  as  in  PDT  II-2  (Figure  B-9),  noncohesive  ejecta 
from  the  surface  surged  up  in  rounded  billows  (but  at  lower  angles  than  in 
PDT  II-2),  moving  outward  with  slight  expansion  (Figures  B-4  and  B-5). 


'k 

Overhead  16-mm  film  made  by  the  Williamson  Aircraft  Company  and  provided 
by  Field  Command,  DNA. 
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FIGURE  B-8  OVERHEAD  VIEW  OF  DICE  THROW  CLOUD  AT  T  +  0.8 


Similar  to  nuclear  "base  surge,"  this  "dirt  surge"  was  caused  by  gust- 
lofted  debris  and  some  missiles  (rock  chunks)  falling  back  onto  the  sur¬ 
face.  After  the  rock  chunks  and  heavier  dirt  clods  rapidly  fell  out, 
these  ejecta  probably  comprised  most  of  the  late-time  dust  cloud  material. 
By  T  +  6  s  (Figure  B-6),  this  ejected  surge  reached  an  altitude  of  only 
100  to  150  m,  much  less  than  the  300  m  that  the  PDT  II-2  spires  had 
attained  by  5  s.5 

The  center  protuberance,  presumably  composed  principally  of  detona¬ 
tion  products,  reached  an  almost  stationary  height  by  about  T  +  6  s,  al¬ 
though  measurements  from  two  views  (R-l  and  PR)  give  evidence  of  a  5-to-7-m 
oscillation  during  the  time  span  T  +  4  s  to  8  s.  This  is  consistent  with 
the  velocity  field  calculations  provided  by  AFWL,  which  now  utilize  a 
newly  developed  "special  rezone  technique"  of  the  HULL  code.  Figure  B-10 
displays  the  HULL  velocity  field  calculation  for  DICE  THROW^  at  T  +  5  s 
superimposed  to  scale  on  the  R-l  view  of  DICE  THROW.  This  computation 
predicts  a  downdraft  in  the  center  of  the  cloud  that  is  part  of  a  "reverse" 
vortex  in  the  cloud.  Thus,  the  initial  upward  movement  of  the  cloud 
center  is  rapidly  stopped;  moreover,  according  to  the  calculation,  a  normal 
vortex  is  beginning  to  form  outside  the  dust-laden  region. 


Cratering 


A  study  of  crater  geometry  is  necessary  in  order  to  determine  the 

approximate  mass  of  lofted  dust  and  debris.  The  apparent  crater  formed 

± 

by  DICE  THROW  was  about  54  m  in  diameter  and  8  m  deep.  Figures  B-4(b)  and 
and  B-8  show  the  scaled  crater  outline  overlayed  on  photographs  taken  from 
the  PR  site  and  from  above  GZ,  respectively. 


Descriptions  of  this  technique  were  published  in  Appendix  C  of  Bollen 
et  al.,s  in  the  form  of  a  memorandum  authored  by  Dr.  Charles  Needham 
of  AFWL  and  in  the  Proceedings  of  the  DICE  THROW  Symposium  as  well. 

^The  HULL  code  calculations  for  DICE  THROW  presented  in  this  report  were 
kindly  provided  by  Major  G.  P.  Ganong  of  AFWL. 

Cratering  data,  which  were  originally  supplied  to  us  by  Field  Command, 
DNA,  have  been  published  in  the  Proceedings  of  the  DICE  THROW  Symposium.1 
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FIGURE  8-10  DICE  THROW  CLOUD  COMPARED  TO  HULL  VELOCITY 
FIELD  CALCULATION  AT  5  s 


As  discussed  earlier,  water  table  conditions  for  the  DICE  THROW 
event  probably  prevented  the  ejection  of  as  much  dense  material  into  the 
air  as  occurred  in  PDT  II- 2.  In  addition,  the  local  geologic  setting  was 
deliberately  modified  in  order  to  reduce  the  number  of  large  particles  or 
"missiles"  for  DICE  THROW,  in  the  hope  of  reducing  ejecta-caused  damage 
to  the  equipment  used  in  the  experiments.  Since  the  ejected  dirt  in  DICE 
THROW  was  rather  noncohesive  gypsum  sand,  one  would  expect  a  smaller 
crater  in  DICE  THROW  than  in  PDT  II-2.  This  should  then  result  in  a  lower 
cratering  efficiency,  measured  by  the  cubic  feet  of  material  permanently 
displaced  from  the  crater  region  per  ton  of  TNT  yield. 

In  order  to  determine  the  maximum  amount  of  crater  ejecta  lofted 
inco  the  air,  it  is  necessary  to  estimate  the  amount  of  soil  missing  from 
the  crater.  This  is  done  by  calculating  the  volume  of  the  apparent 
crater  and  reducing  it  by  the  volume  of  the  dirt  forming  the  lip.  The 
crater  and  lip  profiles  were  measured  shortly  after  the  test.  Figure  B-ll 
shows  the  cumulative  volume  as  a  function  of  the  distance  from  the  crater 
center,  numerically  integrated  from  the  measured  profile. 

From  the  data  presented  in  Figure  B-ll,  it  appears  that  the  "missing" 

3 

volume  of  crater  ejecta  asymptotically  approaches  1560  m  .  This  volume 

3 

is  significantly  less  than  the  cratering  volume  of  4250  m  predicted  for 

3 

DICE  THROW;  that  prediction  was  based  on  an  estimate  of  8.5  m  of  ejecta 
per  ton  of  TNT  and  on  the  500-ton  TNT  equivalence  for  the  DICE  THROW  event. 
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FIGURE  B-11  CUMULATIVE  CRATER  VOLUME  vs  DISTANCE  AWAY  FROM 
CRATER  CENTER 


O 

For  an  average  soil  density  of  approximately  2  g/cm  ,  the  maximum  total 
mass  of  DICE  THROW  ejecta  lofted  an  appreciable  distance  into  the  air  is 
estimated  to  be  3. A  kt.  This  value  is  about  half  of  our  predicted  value 
of  8  to  9  kt.  Much  of  the  discrepancy  may  result  from  reducing  the  ap¬ 
parent  crater  volume  by  the  volume  of  material  in  the  lip.  Although  the 
lip  material  may  have  been  part  of  the  "cloud"  at  early  times,  the  smaller 
value  is  better  at  late  times  because  it  more  fairly  represents  the  amount 
of  soil  carried  away  in  the  cloud.  The  value  chosen  for  the  mass  aloft 
should  be  determined  by  the  time  segment  being  studied. 

5.  Intermediate  Time--T  +  6  s  to  T  +  15  s 

In  this  epoch,  the  darker  protuberance  in  the  center,  produced  by 
the  material  used  in  the  detonation,  disappears  from  view  as  the  surround¬ 
ing  dust  cloud  rises  above  it  (Figures  B-12  and  B-13).  This  behavior  is 
consistent  with  the  predicted  downdraft  in  the  center  of  the  cloud  (see 
Figure  B-10).  Photographs  taken  from  the  aircraft  at  6  and  15  s  [Figures 
B-14(a)  and  (b)]  are  also  consistent  wit!  an  initial  center  downdraft, 
exterior  vortex,  and  low  velocities.  The  HULL  velocity  field  calculation 
at  10  s,  shown  superimposed  on  the  cloud  outline  in  Figure  B-15,  shows 
only  the  normal  vortex  direction--the  close-in  reverse  vortex  has 
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FIGURE  B-15  DICE  THROW  CLOUD  COMPARED  TO  HULL  VELOCITY 
FIELD  CALCULATION  AT  10  s 


dissipated.  But  this  comparison  suggests  that  no  material  was  injected 
into  the  region  where  the  normal  vortex  formation  was  predicted.  The 
general  appearance  of  the  cloud  resembles  a  "doughnut"  with  a  "hollow"  in 
the  center  by  15  s  (Figure  B- 14(b)]. 

By  about  T  +  10  s,  it  is  readily  evident  that  the  cloud  had  begun  to 
move  appreciably  to  the  northwest  or  away  from  the  phase  repeater,  to  the 
left  of  R-l,  and  to  the  right  of  the  line  of  sight  from  North  Oscuro  Peak. 
Before  T  +  15  s,  the  residual  dust  cloud  maintained  an  overall  shape  simi¬ 
lar  to  that  at  early  times  except  for  the  absence  of  the  upward-protruding 
structure  of  detonation  material.  During  this  time  sequence  the  main 
part  of  the  cloud  steadily  rose  and  expanded  around  the  dome  of  detonation 
products  in  the  center.  At  about  T  +  15  s,  particularly  as  viewed  from 
the  R-l  and  NOP  sites,  a  peaestal-like  "stem"  started  to  form  as  the 
lower  part  of  the  cloud  above  the  ground  dust  layer  became  "pinched-of f" 
[Figure  B-16(a)  and  (c)]. 

6.  Late  Time--After  T  -F  15  s 

If  it  had  behaved  in  a  way  similar  to  a  nuclear  or  tangent-sphere 
HE  detonation  after  T  +  15  s,  the  DICE  THROW  cloud  should  have  been  af¬ 
fected  by  cool  air  rushing  in  under  the  rising  vortex  ring  and  bringing 
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in  residual  dust  or  "sweep-up”  material.  (Figures  B-17,  B-18,  and  B-19 
show  such  characteristics  for  nuclear  and  tangent-sphere  TNT  detonations 
after  T  +  15  s.)  However,  according  to  HULL  calculations,  which  the 
photographic  data  support,  the  DICE  THROW  vortex  was  outside  the  center  of 
the  dust  cloud  and  was  quite  weak.  Figure  B-20  shows  both  the  HULL  linear- 
scale  (superimposed  on  ♦•he  cloud  outline)  and  the  log-scale  velocity  field 
calculations  at  T  +  20  s.  The  maximum  predicted  velocity,  which  was 
5.6  m/s,  steadily  decreased  with  time  after  detonation.  The  log-scale 
velocity  field  illustrates  the  extent  of  the  circular  motion  of  the  com¬ 
puted  vortex.  In  addition,  there  was  no  longer  a  predicted  downdraft  at 
the  center  of  the  cloud. 

A  stem  eventually  developed  (Figures  B-21  and  B-22)  that  was  bulky, 
compared  to  tangent-sphere  HE  and  nuclear  detonations.  By  T  +  45  s,  a 
very  evident  updraft  began  to  influence  the  lofting  of  sweep-up  and  crater 
ejecta  material  up  the  "stem"  (Figure  B-23).  In  a  nuclear  cloud  case,  the 
dust  traveling  up  the  stem  becomes  entrained  in  the  rising  vortex  ring 
forming  the  traditional  "mushroom"  cap;  however,  for  DICE  THROW,  the  dust, 
which  appeared  to  move  according  to  AFWL  calculations,  did  not  seem  to 
become  entrained  in  a  vortex.  Figure  B-24  compares  velocity  field  calcu¬ 
lations  with  the  DICE  THROW  cloud  at  T  +  30  s  and  T  +  45  s.  The  maximum 
(updraft)  velocity  slowly  increased,  sweeping  the  dust  up,  but  the  pre¬ 
dicted  weak  vortex  ring  remained  outside  of  the  visible  cloud.  We  believe 
that  entrainment  of  dust  in  the  vortex  did  not  occur.  Rather,  the  late¬ 
time  morphology  for  DICE  THROW  is  best  described  as  being  an  amorphous 
and  slowly  rising  stem,  composed  of  several  convective  cells  displaying 
their  own  circulation  patterns,  rather  than  a  rapidly  rising  vortex  ring 
with  dust  entrainment. 

PDT  II-2  did  not  exhibit  such  classical  "mushroom  shape"  properties 
as  strong  sweep-up  winds,  seem,  or  vortex  ring  either  (Figure  B-25). 
Therefore,  because  of  its  same  type  and  shape  of  charge,  DICE  THROW  was 
not  expected  to  develop  such  features.  It  was  anticipated  that  DICE  THROW 
would  possess  cloud  morphology  similar  to  that  from  buried  (as  deep  as 
17  ft)  nuclear  detonations--e.g. ,  Buster  Jangle:  Uncle  (Figure  B-26).5 
On  the  contrary,  DICE  THROW  appeared  less  amorphous  than  PDT  11-2  and 
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FIGURE  B-17  NUClEAP.  DETONATIONS  AT  T 


(e)  T  40  * 

FIGURE  B-18  100-T  TNT,  PRE-DICE  THROW  II-1 
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FIGURE  B-21  DICE  THROW  AT  20  s 
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FIGURE  B-23  DICE  THROW  AT  45  s 
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FIGURE  B-24 


DICE  THROW  CLOUD  COMPARED  TO  HULL  VELOCITY  FIELD 
CALCULATIONS 


FIGURE  B-26  1.2-kt  NUCLEAR  AT  -17  ft,  BUSTER  JANGLE:  UNCLE  1951,  T  +  20  s 


seemed  to  exhibit  some  mushroom  cloud  and  stem-type  features.  To  a  more- 
than- anticipated  degree,  it  did  resemble  surface  or  shallow  subsurface 
nuclear  detonations  such  as  Johnny  Boy  or  Buster  Jangle:  Sugar,  and,  since 
a  stem  formed,  it  had  more  similarity  to  TNT  tangent  sphere  detonations 
such  as  PDT  II-l  (Figure  B-18)  and  Mixed  Company  III  (Figure  B-19)  at 
times  after  T  +  20  s  than  expected.  However,  the  most  notable  difference 
between  PDT  II- 2  and  DICE  THROW  was  the  lack  of  heavy  ejecta  spires  in 
the  latter. 

A  strong  wind  shear  rapidly  deformed  the  DICE  THROW  cloud.  The  lower 
parts  of  the  cloud  moved  to  the  northwest,  but  as  the  cloud  top  rose,  it 
was  swept  around  to  the  southeast.  This  distortion  is  readily  apparent 
in  the  photographs  for  T  +  1  m  and  T  +  1  m  30  s  (Figure  B-27),  and  makes 

it  difficult  at  late  times  to  compare  the  AFWL  flow  fields  with  the  actual 

* 

cloud  shape. 

Figure  B-28  compares  the  cloud  at  one  minute  to  the  AFWL  velocity 
field  calculation.  The  predicted  vortex  remains  outside  of  the  visible 
cloud  and  is  still  very  weak.  Although  the  predicted  updraft  velocity  in 
the  cloud  center  has  increased  to  9.5  m/s,  it  is  nevertheless  substantially 
lower  than  that  associated  with  tangent  sphere  HE  detonations.  Consistent 
with  the  predicted  weak  vortex  velocities,  the  material  lofted  appeared 
to  remain  as  a  stem. 

As  seen  from  the  PR  site  [Figure  B-29(b)],  the  main  cloud  began  to 
break  away  from  its  stem  at  around  T  +  2  m.  This  separation  was  well 
developed  by  T  +  4  m,  by  which  time  dissipation  of  both  the  cloud  and 
stem  has  set  in  (Figure  B-30).  From  personal  observations,  after  T  + 

4  min,  the  main  part  of  the  cloud  eventually  moved  to  the  southeast, 
rising  and  gradually  dispersing  up  and  over  the  South  Oscuro  ridge.  The 


Ganong  et  al.  ,11  were  quite  successful,  however,  in  reproducing  the 
appearance  of  the  DICE  THROW  cloud  at  late  times  by  modifying  tracer 
particle  positions  in  their  HULL  calculations  according  to  the  wind 
profile. 
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FIGURE  B-30  DICE  THROW  AT  T  +  4  min 


lower  portions  of  the  cloud  continued  to  be  blown  to  the  northwest,  and 
a  large  part  of  the  Jornada  de  Muerto  basin  took  on  a  hazy  appearance. 

7.  DICE  THROW  Dimensional  Measurements 

Scaled  measurements  of  DICE  THRCW  morphology  were  undertaken  to  sup¬ 
port  interpretation  of  the  rf  data.  These  measurements  were  made  by 
measuring  directly  from  the  film  records  (or  magnified  prints  thereof) 
and  then  factoring  in  the  appropriate  variables  such  as  distance  between 
the  camera  and  cloud  at  Tq  and  focal  length  of  the  lens.  In  some  cases, 
where  distance  was  questionable,  known  distances  between  objects  in  the 
field  of  view  close  to  ground  zero  were  utilized  to  determine  the  mag¬ 
nification  factor.  Accuracy  is  estimated  to  be  with  ±5  m;  this  is 
established  through  comparison  of  specific  dimensions  measured  from  dif¬ 
ferent  film  records  at  the  same  time  after  detonation,  taking  into  account 
variables  in  the  aspect  view  of  cloud  movement  and  morphology.  An  expo¬ 
sure  of  an  Accutron  clock  face  was  included  in  each  picture  frame  for 
timing. 

Figures  B-31  through  B-33  give  the  measured  DICE  THRCW  dimensions  of 
cloud  height,  cloud  width,  and  stem  width  as  a  function  of  time  after 
detonation  taken  from  photographic  data,  recorded  in  cameras  at  the  R-l, 
PR,  and  NOT  sites,  respectively.  In  each  case,  the  cloud  height  was 
measured  from  the  ground  to  the  highest  visible  point  on  the  cloud,  width 
was  determined  to  be  the  widest  extent,  and  stem  width  was  measured 
across  its  vertical  center.  In  addition,  the  slope  of  the  rate  of  cloud 
rise  is  identified  and  indicated.  Of  particular  interest  are  the  cloud 
height  measurements.  By  comparing  the  different  views  with  each  other 
and  with  the  photographs  in  the  previous  section,  one  can  form  a  general 
idea  of  gross  cloud  movement  caused  by  prevailing  winds. 

Up  to  about  T  +  15  s  the  cloud  height,  viewed  from  the  three  sites, 
was  fairly  even.  Thereafter  the  height  measured  from  the  R-l  photographic 
record  (Figure  B-31)  attained  higher  levels  up  to  about  T  +  1  min.  Data 
measured  from  the  PR  (Figure  B-32)  and  NOP  (Figure  B-33)  sites  were 
fairly  even  at  T  +  20  s;  thereafter,  the  cloud  height:  from  the  NOP  data 
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FIGURE  B-31  MEASURED  CLOUD  HEIGHT,  WIDTH,  AND  STEM  WIDTH  vs  TIME  FROM  RECEIVER  1  CAMERA  SITE 
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FIGURE  B-33  MEASURED  CLOUD  HEIGHT,  WIDTH,  AND  STEM  WIDTH  vs  TIME  FROM  NORTH  OSCURO  PEAK  CAMERA 


increased  more  than  that  from  the.  PR  data  and  rose  above  that  from  the 
R-l  data  at  T  +  90  s.  Rise  rates  for  both  R-l  and  NOP  site  data  were 
fairly  even  in  the  10-to-90-s  regime  (13.1  to  13.2  m/s)  and  higher  than 
the  data  exhibited  by  the  PR  records  (11.6  m/s).  The  discrepancies  were 
most  likely  due  to  low-altitude  prevailing  winds  carrying  the  cloud  away 
from  the  PR  site  and  somewhat  closer  to  the  R-l  site  to  the  northwest. 

Of  course,  as  the  cloud  rose,  the  top  edge  moved  relatively  closer  to 
the  top  of  North  Oscuro  Peak,  thus  explaining  the  greater  height  recorded 
from  the  NOP  site  at  T  +  90  s  than  from  the  R-l  site. 

The  measured  cloud  expansion,  as  seen  from  the  PR  and  NOP  sites,  was 
somewhat  faster  than  from  the  R-l  site.  Later,  after  T  +  30  s,  data  from 
NOP  showed  the  greatest  expansion,  while  that  from  the  PR  site  revealed  the 
slowest.  The  prevailing  winds  that  moved  the  DICE  THROW  cloud  away  from 
the  PR  site  can  explain  part  of  this  discrepancy.  However,  a  more  import¬ 
ant  effect  is  that  of  the  change  in  the  appearance  of  the  cloud,  when 
viewed  from  different  directions.  Figure  B-34  gives  height  and  width 
dimensions  for  DICE  THROW  averaged  from  the  R-l,  PR,  and  NOP  sites;  the 

average  rise  rate  between  T  +  10  and  T  +  120  s  was  13  m/s.  Figure  B-35 

* 

presents  the  corresponding  PDT  II-2  quantities.  Between  T  +  10  and  T  + 

60  s,  the  PDT  II-2  cl  ud  had  a  rise  rate  of  only  7.6  m/s  compared  to 
13  m/s  for  DICE  THROW  during  the  same  time  period,  and  by  T  +  60  s  the 
PDT  II-2  cloud  had  attained  a  height  of  only  600  m  compared  to  about  800  m 
for  DICE  THROW.  The  ratios  of  rise  rate  and  altitude  attained  (1.7  and 
1.3,  respectively)  are  consistent  with  a  cube-root-of -yield  scaling  rule 
[(628/120)l/3  =  1.7]. 

The  effect  that  prevailing  winds  had  on  displacing  the  DICE  THROW 
cloud  away  from  ground  zero  as  seen  from  the  R-l  site  is  graphically 
presented  in  Figure  B-36,  Up  to  about  1200  m  of  height  (90  s  after 
detonation),  the  cloud  drifted  relatively  uniformly  to  the  northwest. 

At  higher  altitudes,  however,  winds  moving  more  to  the  south  or  southeast 

•k 

Measurements  were  made  from  photos  provided  by  A.  D.  Thornbrough, 

Sandia  Labs. 


FIGURE  B  35  MEASURED  PDT  II-2  CLOUD  HEIGHT 
AND  WIDTH  vs  TIME 

pushed  the  north  side  of  the  cloud  top  along  with  them.  The  center  and 
opposite  side  of  the  cloud  were  more  gradually  affected  by  this  change 
in  wind  vector  because  of  the  previous  displacement  in  the  opposite 
direction  and  the  expansion  of  the  cloud.  Thus,  as  noted  in  the  photo¬ 
graphs,  the  lower  portions  and  stem  of  the  DICE  THROW  cloud  maintained  a 
northwest ?rly  movement  while  the  upper,  main  part  of  the  cloud  was  blown 
to  the  southeast,  gradually  rising  and  dispersing  up  over  the  South 
Oscuro  ridge. 

Figure  B-37  displays  the  DICE  THROW  cloud  volume  as  a  function  of 
time  as  measured  from  the  film  records  of  the  three  major  camera  stations, 
R-lj  PR>  and  NOP,  The  volume  data  from  the  PR  site  were  corrected  slightly 
(increased)  to  allow  for  the  increase  in  distanct  with  time  between  the 
camera  and  the  cloud.  These  data  can  be  used  to  estimate  the  average 
dust  density  in  the  cloud  as  a  function  of  time.  Figure  B-38  shows  the 
resulting  estimated  average  dust  density  between  T  +  1  and  T  +  100  s  for 
a  constant  mass  aloft  of  3.4  kt,  the  value  obtained  from  crater  measure¬ 
ments  described  above.  Also  shown  is  the  average  dust  density  for  T  +  1 
to  i  +  6  s  for  10  kt  aloft.  During  the  initial  phase  when  most  of  the 
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CLOUD  DISPLACEMENT  — m 

FIGURE  8-36  DICE  THROW  CLOUD  EDGE  DISPLACEMENT  FROM  GZ  (R-1  Site) 


rapid  fallout  probably  occurred,  the  average  dust  density  in  the  DICE 

A  —  3  3 

THROW  cloud  was  in  the  range  10’  to  10  g/cm.  According  to  this 

-5  3 

model,  the  average  density  had  dropped  to  about  10  g/cm  by  T  +  100  s. 

8.  HULL-Computation  and  DICE  THROW  Dust  Cloud  Predictions 

Figures  B-39  and  B-40,  which  are  reproduced  from  Ganong  et  al.,11 
summarize  some  post-event  HULL-code  calculations.  Figure  B-39  presents 
the  mass  aloft  as  a  function  of  time  for  two  particle  size  distributions. 
Figure  B-40  shows  the  corresponding  values  of  Y,  which  was  defined  in 
Eq.  (A-7).  Note  that  this  DICE  THROW  distribution  is  not  the  same  as  the 
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FIGURE  B-38  AVERAGE  DUST  DENSITY  IN  THE  DICE  THROW  CLOUD  INFERRED 
FROM  VOLUME  MEASUREMENTS 


one  approximated  by  Eq.  t.A-14),  which  has  features  of  both  the  DICE  THROW 
and  hard  rock  distributions.  These  results  we  obtained  by  assuming  that 
10  kt  of  ejecta  were  initially  deposited  in  the  volume  occupied  by  the 
actual  DICE  THROW  cloud  at  T  +  1  s.  The  mass-aloft  calculation  agrees 
fairly  well  with  the  cratering  volume  and  "missing"  mass  results  of 
Section  4  of  this  appendix. 
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FIGURE  B-39  MASS  OF  DUST  ALOFT  vs  TIME  FOR  DICE 
THROW  AND  HARD  ROCK  PARTICLE  SIZE 
DISTRIBUTIONS 


9.  Comparison  Between  Predicted  and  Actual 
Propagation-Path  Cloud  Penetration  Points 

One  of  the  important  sources  of  information  for  the  criteria  used  to 
determine  the  physical  layout  of  this  experiment  was  the  PDT  II-2  cloud 
development.  The  PDT  II-2  event  was  a  120-ton  ANFO  detonation  that  had 
exactly  the  same  spherically  capped  cylinder  configuration  as  the  DICE 
THROW  Main  Event  (PDT  II-2  was  the  culmination  of  the  DICE  THRCW 
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FIGURE  B-40  ANFO  DUST  CLOUD  CASE  BY  PARTICLE  SIZE  DISTRIBUTION  POWER  (7)  vs 
TIME  FOR  DICE  THROW  AND  HARD  ROCK  PARTICLE  SIZE  DISTRIBUTIONS 


configuration  development).  Because  the  PDT  II-2  dust  cloud  evolved  in 
a  much  different  way  than  that  of  the  PDT  II-l  event,  which  used  a  100- 
ton  TNT  tangent  sphere  configuration,  PDT  II-2  was  believed  to  be  the 
best  guide  for  DICE  THROW  planning. 

Figure  B-41  shows  the  six  DICE  THROW  propagation-path  penetration 

poincs  at  several  times  superimposed  on  photographs  of  the  PDT  II- 2  dust 

cloud.  Eoth  the  times  and  dimensions  have  been  scaled  according  to  a 

3 

cube-root-of -yield  rule  (i.e.,  by  */5  ).  If  this  procedure  were  accurate, 
a  series  of  long-duration  occultations  would  have  been  obtained.  However, 
as  seen  in  Figure  B-42,  where  the  actual  DICE  THROW  penetration  points 
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FIGURE  B-41  YIELD-SCALED  PHOTOGRAPHS  OF  PDT  II-2  SHOWING  ANTICIPATED  INTERSECTION  OF  PROPAGATION  PATHS  WITH 
DICE  THROW  PHENOMENA.  Dimensions  and  time  scaled  as  3%/IT  (Photos  courtesy  of  A.  D.  Thornbrough,  Sandia  Labs.) 


FIGURE  B-42  PHOTOGRAPHS  OF  DICE  THROW  SHOWING  ACTUAL  INTERSECTION  OF  PROPAGATION  PATHS 


are  shown  at  about  the  same  times  as  the  scaled  ones  of  Figure  B-41,  the 
predictions  were  quite  wrong.  Two  major  reasons  for  the  disagreement 
were  the  high  winds  during  DICE  THROW  and  the  different  soil  conditions 
between  the  two  sites  (the  PDT  series  took  place  at  the  Queen  15  site 
about  50  km  south  of  the  Giant  Patriot  location). 

Unlike  the  DICE  THRCW  soil,  the  PDT  soil  was  very  wet.  Large  craters 
were  formed  there,  and  numerous  large  (tens  of  centimeters  across)  clay 
"missiles"  were  ejected  and  fell  over  a  large  area.  The  rapidly  develop¬ 
ing  ejecta  spires  that  can  be  seen  in  the  first  three  frames  of  Figure  B-41 
were  apparently  composed  of  those  large  particles.  Concern  that  the  "bombs" 
might  damage  the  equipment  for  various  experiments  arranged  around  the 
DICE  THRCW  detonation  led  to  the  effort  to  suppress  them  by  excavating  the 
area  around  ground  zero,  and  backfilling  it  with  sand.  It  should  not  be 
surprising,  therefore,  that  the  spires  were  not  seen  during  DICE  THRCW. 

The  PDT  II- 2  and  DICE  THROW  clouds  appear  to  have  developed  in  about  the 
same  way,  otherwise,  if  the  deformation  caused  by  wind  shear  is  taken 
into  account.  However,  the  presence  of  many  very  large  particles  at 
early  times  during  PDT  II- 2  would  have  undoubtedly  led  to  disproportionally 
strong  rf  effects. 

10.  Summary 

From  photographic  studies,  the  overall  cloud  morphologies  of  DICE 
THROW  and  PDT  II- 2  do  not  resemble  each  other  enough  to  unequivocally 
state  that  geometric  predictions  for  the  DICE  THRCW  behavior  were  valid. 

The  anticipated  early  development  of  an  intense  cloud  of  ejecta  (spires) 
to  5  s  or  10  s  did  not  characterize  DICE  THROW.  Perhaps  because  of  its 
larger  yield  and  stronger  flow  field,  the  DICE  THROW  main  event  produced 
a  less  amorphous  cloud  than  did  PDT  II-2.  The  general  development  of  the 
DICE  THROW  cloud  and  the  mass  aloft  seems  to  be  consistent  with  the  HULL 
predictions.  It  appears  that  no  dust  was  injected  by  the  detonation  into 
the  region  where  the  normal-direction  vortex  formation  is  predicted.  If 
the  calculations  are  accurate,  no  dust  entrainment  in  the  vortex  region 
ever  took  place,  and  the  visible  cloud  remained  in  a  region  of  mostly 
vertical  flow  and  thus  corresponded  to  a  "stem."  Because  part  of  the 
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reason  for  performing  measurements  during  these  tests  is  to  verify  or 
improve  the  codes,  the  possibility  of  marking  the  region  of  (currently 
visible)  vortex  formation  with  smoke  trails  or  the  like  should  be 
considered. 
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ATTN:  J.  Krlstof 

Santa  Fe  Corp. 

ATTN:  E.  Ortlleb 

Science  Applications,  Inc. 

ATTN:  L.  Llnson 
ATTN:  D.  Hamlin 
ATTN:  D.  Sachs 
ATTN:  E.  Straker 
ATTN:  J.  rtcDougall 
ATTN:  C.  Smith 

Science  Applications,  Inc. 

ATTN:  D.  Dlvls 

SRI  International 


ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


C.  Rlno 
G.  Smith 
G.  Price 

R.  Leadabratid 
W.  Jaye 
W.  Chesnut 
M.  Baron 
A.  Burns 

D.  Nell  son 

R.  Livingston 


Teledyne  Brown  Engineering 
ATTN:  R.  Dellberls 

Tri-Corn,  Inc. 

ATTN:  D.  Murray 

TRW  Defense  &  Space  Sys.  Group 
ATTN:  S.  Altschuler 
ATTN:  R.  Plebuch. 

ATTN:  D.  Dee 

Utah  State  University 
ATTN:  L.  Jensen 
ATTN:  K.  Baker 

Vtsldyne,  Inc. 

ATTN:  J.  Carpenter 
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